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ABSTRACT 

I derive the physical properties of thirty transiting extrasolar planetary systems using a homo- 
geneous analysis of published data. The light curves are modelled with the JKTEBOP code, 
with special attention paid to the treatment of limb darkening, orbital eccentricity, and error 
analysis. The light from some systems is contaminated by faint nearby stars, which if ignored 
will systematically bias the results. I show that it is not realistically possible to account for 
this using only transit light curves: light curve solutions must be constrained by measure- 
ments of the amount of contaminating light. A contamination of 5% is enough to make the 
measurement of a planetary radius 2% too low. 

The physical properties of the thirty transiting systems are obtained by interpolating in 
tabulated predictions from theoretical stellar models to find the best match to the light curve 
parameters and the measured stellar velocity amplitude, temperature and metal abundance. 
Statistical errors are propagated by a perturbation analysis which constructs complete error 
budgets for each output parameter These error budgets are used to compile a list of systems 
which would benefit from additional photometric or spectroscopic measurements. 

The systematic errors arising from the inclusion of stellar models are assessed by using 
five independent sets of theoretical predictions for low-mass stars. This model dependence 
sets a lower limit on the accuracy of measurements of the physical properties of the systems, 
ranging from 1% for the stellar mass to 0.6% for the mass of the planet and 0.3% for other 
quantities. The stellar density and the planetary surface gravity and equilibrium temperature 
are not affected by this model dependence. An external test on these systematic errors is 
performed by comparing the two discovery papers of the WASP-1 1 /HAT-P-10 system: these 
two studies differ in their assessment of the ratio of the radii of the components and the 
effective temperature of the star. 

I find that the correlations of planetary surface gravity and mass with orbital period have 
significance levels of only 3. la and 2.3(7, respectively. The significance of the latter has not 
increased with the addition of new data since Paper II. The division of planets into two classes 
based on Safronov number is increasingly blurred. Most of the objects studied here would 
benefit from improved photometric and spectroscopic observations, as well as improvements 
in our understanding of low-mass stars and their effective temperature scale. 

Key words: stars: planetary systems — stars: binaries: eclipsing — stars: binaries: spectro- 
scopic — stars: fundamental parameters 



1 INTRODUCTION 

The study of extrasolar plane ts is scientifically and c ulturally im- 
portant, and after a late start ( iMavor & Ouelozll 19951) the number 
of known planets is escalating rapidl)^. Transiting planets are the 
crown jewels of this population as, with the exception of our own 
Solar system, they are the only planets whose masses and radii are 
directly measurable. In addition to this, it is possible to put con- 
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straints on the properties of their atmospheres, in which much in- 
teresting physics occurs, through measurements of the depths of 
transits and occultations at different wavelengths. 

Whilst transiting extrasolar planets (TEPs) offer unique sci- 
entific possibilities, their study involves several complications. The 
most significant is that it is not in general possible to measure the 
mass and radius of a planet through basic observations alone. Ad- 
ditional constraints are needed, and are usually provided by forcing 
the properties of the host stars to match theoretical expectation^- 

^ In principle, astrometric observations could either replace radial velocity 
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This introduces not only a model dependence (i.e. systematic er- 
ror), but also the possibility of inconsistent results if different the- 
oretical predictions are used for some TEPs. Systematic errors can 
blur any distinctions between planets, making it hard to pick out 
discrete groups of TEPs from the varied general population. 

This systematic error cannot be abolished, but it can at least be 
standardised. In this series of papers I am analysing the known tran- 
siting systems using rigorously homogeneous methods, with the 
aim of removing the systematic differences in measurements of the 
physical properties of TEPs. The resulting physical properties are 
therefore statistically compatible, and any structure in distributions 
of parameters is maximised. 

In Paper I l lSouthworthll2008l) I analysed the light curves of 
the fourteen transiting systems for which high-precision photome- 
try was then available, paying careful attention to the role of limb 
darke ning and to the es timation of comprehensive errorbars. Pa- 
per II [SouthworthllioO^ used these results plus the predictions of 
three different theoretical stellar models to measure the physical 
properties of the fourteen TEPs. In this work I broaden the analysis 
to thirty TEPs and five sets of theoretical stellar models, resulting 
in improved statistics and better systematic error estimates. 

There are a few homogeneous analyses of transiting systems 
available in the literature. A good analysis of 23 systems was pre- 
sented bv [Torres et al.l ( |2008| , hereafter TWH08), but these authors 
tried only two different theoretical model sets and did not assign 
systematic errors to their results. Analogously, such work would 
also benefit from homogeneous analysis of the spectra of the host 
stars in order to put their effective temperature and chemical abun- 
dance measuremen ts on a consistent s cale. Steps towards this goal 
were p ioneered by Valenti & Fischer ( 2005) a nd are being co ntin- 
ued bv lAmmler-von Eiff et al. (2009) and Ghezzi et alj ( I2OI0I) , but 
a homogeneous study of the host stars of all known TEPs is not 
currently available. 

In Sect.[2]l present the methods used to analyse the light curves 
of the thirty TEPs included in this work. Sect.[3] discusses the five 
theoretical stellar model sets and their application to determining 
the physical properties of the TEPs. Sect.|4]presents the new results 
for these objects, Sect.|5] discusses the influence of systematic er- 
rors due to the use of theoretical models, and in Sect.[6]l summarise 
the physical properties of the known TEPs and explore correlations 
between various parameters. Those readers interested in the gen- 
eral properties of TEPs rather than specific systems can skip Sect.|4] 
without problem. 



2 LIGHT CURVE ANALYSIS: JKTEBOP 

I have mod elled the light curves o f each TEP using the JKTE- 
8010 code jSouthworth et alj[2004al lbh. JKTEBOP grew out of the 
original E BOP program written for eclipsing binary star systems 
(Popper & Etzel 1981; Etzel 1981) and implementing the NDE 
model (iNelson & Davis 1972.). JKTEBOP uses biaxial spheroids to 
model the component stars (or star and planet) so allows for depar- 
tures from sphericity. The shapes of the components are governed 
by the mass ratio, q, although the results in this work are all ex- 
tremely insensitive to the value of this parameter. 



The main parameters of a JKTEBOP fit are the orbital inclina- 
tion, j, and the fractional radii of the two starfl ta and rb. The 
fractional radii are defined as 



rA = 



Ra 



Rh 

a 



(1) 



where Ra and i?b are the stellar and planetary radii and a is the 
orbital semimajor axis, t-a and rb correspond to radii of spheres 
of the same volume as the biaxial spheroids. In JKTEBOP the frac- 
tional radii are re-parameterised as their sum and ratio: 
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because these are only weakly correlated with each other. In general 
the orbital period, Porb, is taken from the literature and the time 
of transit midpoint, To, is included as a fitted parameter in each 
JKTEBOP run. 



2.1 Treatment of limb darkening 

The limb darkening (LD) of the star is an important 'nuisance pa- 
rameter' affecting transit light curves which can be parametrised 
using any of five LD laws in JKTEBOP. Wherever possible the LD 
coefficients are included as fitted parameters, but when there is in- 
sufficient information for this the coefficients are fixed at theoreti- 
cal values. For each light curve I have obtained solutions with each 
of the five LD laws (see Paper I for their definition) and with both 
LD coefficients fixed, with the linear coefficient fitted and the non- 
linear coefficient fixed (hereafter referred to as 'LD fit/fix'), and 

with both coefficients fitted. 

T heore tical LD coefficients have been taken from'Van Hamrn^ 

( Il993l) .lciaret (2000, 2004a) and Claret & Hauschildt (2003). The 
tabulated values have been bilinearly interpolated, using the JKTLD 
codfl to the known effective temperature (Tcft) and surface gravity 
(log g) of the star. I find that there is usually a spread of 0.1-0.2 in 
the theoretical LD coefficients for cool stars, so when the nonlinear 
LD coefficient is not included as a fitted parameter it is perturbed 
in the Monte Carlo simulations by ±0.1 on a flat distribution to 
account for this. The dependence on theoretical calculations in this 
case is still exceptionally small, because the linear and nonlinear 
coefficients of the LD laws are highly correlated with each other 
( Southworth et al. 2007a). 

Once solutions have been obtained for the five LD laws, the fi- 
nal result is calculated by taking the weighted means of the param- 
eter values for the four two-parameter LD laws (i.e. the linear law 
is not used). The parameter errorbars are taken to be the largest of 
the individual errorbars (see below) plus a contribution to account 
for scatter of the parameter values from different LD law solutions. 



2.2 Error analysis 

For each solution I run 1000 Monte Carlo simulations 
dSouthworth et al.ll2004d,l2005h to provide robust estimates of the 
la statistical errorbars. The starting parameter values are perturbed 



measurements, or augment them and thus provide the missing constraint, 
but this has not yet been achieved in practise. 

JKTEBOP is written in FORTRAN77 and the source code is available at 
|http : / / www .astro ■ keele ■ ac ■ uk/~jkt /codes / jktebop ■ html| 



* Throughout this work stellar parameters are indicated by a subscripted 
'A' and planet parameters by a subscripted 'b', to conform to lAU nomen- 
clature. 

^ JKTLD is written in FORTRAN77 and the source code is available at 
|http : //www .astro ■ keele ■ ac ■ uk/^ jkt/ codes / jkt Id . html 
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for each simulation to avoid sticking artificially close to the origi- 
nal best fit. If the reduced of the fit, x^, is greater than unity, the 
Monte Carlo errorbars are multiplied by ^J'xv to account for this. 

Monte Carlo simulations do not fully account for the presence 
of correlated ('red') noise, which is an unavoidable reality in high- 
precision light curves of bright stars. I theref ore also run a resid- 
ual permutation (or "prayer bead") algorithm ^ Jenkins et al.ll20o3) 
with the quadratic LD law. If there is significant correlated noise 
the residual-permutation errorbars will exceed the Monte-Carlo er- 
rorbars. I then take the larger of the two error estimates to represent 
the final errorbars of the photometric parameters. 



2.3 Orbital eccentricity 

Some TEPs have a non-circular orbit which must be accounted for 
in the light curve analysis. Orbital eccentri city is very dif ficult to 
detect from the shape of a transit light curve jKippindioog) but can 
have a significant effect on the resulting parameters (for an example 
see Sect. |4. 131 . Non-circular orbits normally become apparent from 
radial velocity (RV) measurements of the parent stars. These RVs 
can then be used to determine the eccentricity, e, and the longitude 
of periastron, uj, of the binary orbit, 

JKTEBOP has been modified to account for orbital eccentricity 
by including the possibility of specifying values for either e and 
u) or the combinations ecoso; and esino;. These values and their 
uncertainties are then simply treated as extra observations, and e 
and ui (or their combinations) are included as fitted parameters. In 
this way the uncertainties in e and ui are correctly propagated into 
the errorbars in the other photometric parameters. I prefer to work 
with e cos u) and e sin uj rather than e and u, because t he latter two 
quant ities are strongly correlated with each other (e.g. lBruntt et al.l 
l2006l) . 
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Figure 1. Plot of the vaiiation of pai'ameters fitted to a set of synthetic 
transit light curves with 1 mmag of Gaussian noise added. The synthetic 
datasets were generated for a range of third light values but fitted with the 
assumption of L3 = 0. Dotted lines show the input parameter values for 
the synthetic light curve calculations. 



2.4 Contaminating ligtit 

It is possible for additional light to contaminate photometric obser- 
vations of transiting planets. Any extra light from nearby faint stars 
will dilute the transit depth, causing a systematic error in the light 
curve parameters. This idea is becoming more important for several 
reasons. Firstly, the CoRoT satellite has a large point spread func- 
tion (PSF) which usually contains a number of stars aside from the 
one hosting a transiting planet. Secondly, observations using tele- 
scope defocussing are potentially more susceptible to contaminat- 
ing light. Thirdly, iDaemgen et alj 1^09) have detected faint com- 
panions to three transiting systems (TrES-2, TrES-4 and WASP-2) 
from ground-based high-resolution observations of fourteen TEPs 
obtained with a lucky-imaging camera. These companions could be 
bound to their respective transiting systems, or may just be aster- 
isms. 

Temporarily ignoring the orbital ephemeris (Porb and To), 
there are three main observables in a transit shape: its depth, dura- 
tion, and the duration of totality (Paper I). From transit light curves 
we measure three quantities, which in the case of JKTEBOP are ta, 
rb and i. It is therefore expected to be impossible to fit directly for 
contaminating light, as this would require measuring four indepen- 
dent parameters using only three observables. This expectation will 
now be verified. 



2.4.1 The effect of third light 

I have explored the possibility of measuring contaminating light 
by simulating a set of light curves with reasonable parameters 
(r-A + Ti, ~ 0.1, k — 0.1, i = 87°). I added contaminating light 
(by convention referred to as 'third light' and expressed as a frac- 
tion of the total system light) by amounts ranging from L3 = 
to 0.75 in steps of 0.05. These were transformed into typical good 
ground-based light curves by retaining approximately 400 points 
within each transit and adding a Gaussian scatter with standard de- 
viation 1 mmag. These synthetic light curves were then fitted with 
JKTEBOP under the assumption that La = 0. The results (Fig.O 
show that the presence of is results in systematic overestimates 
of TA and underestimates of rb and i. The bottom panel of Fig.[T] 
shows that the quality of the fit does not get worse as L3 increases. 
Unaccounted third light therefore biases the resulting parameters 
without being detectable through its impact on the quality of the 
fit. 

As a second test I modelled the same synthetic light curves 
again, this time fitting for third light. The resulting values of L3 
are shown in Fig.|2]and are extremely scattered as well as biased to 
smaller values. JKTEBOP deliberately does not restrict photometric 
parameters to physically realistic values (e.g. L3 ^ 0), to avoid 
statistical biases in the uncertainties arising from Monte Carlo sim- 
ulations. Fig.|2] demonstrates that there is a very small amount of 
information on L3 in a good ground-based light curve, but that this 
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Figure 2. Plot of fitted versus input values of third light for the same light 
curves as in Fig.[T] but with L3 included as a fitted parameter. The dotted 
line shows parity. Note the large scale on the j/-axis. 



information is far too sparse to be useful. Fig.[3]shows the resulting 
values of the other main parameters: rb and i are biased towards 
lower values and there is no trend visible in the sizes of the residu- 
als. 

In reality a large value of L3 is not expected because such 
a bright star would show up in the spectroscopic observations 
of a transiting system. Fainter stars can be found via high- 
reso lution ima ging if they are slightly away from the planet host 
star jPaemgen et al. 2009). But it may never be possible to rule out 
the presence of a much fainter star (L3 < 5% depending on the 
quality of the spectroscopic observations) which almost exactly co- 
incide with the planet host. As a guide, 5% of third light can be 
compensated for by increasing va by 1%, and decreasing rb by 
2% and i by 0. 1 ° . It can therefore change the derived radius of the 
planet by several percent. 




0.2 0,4 

hput third light value 

Figure 3. Plot of the variation of parameters fitted to the same light curves 
as in Fig.[T] but in this case with L3 included as a fitted parameter. The 
dotted lines show parity. 



are within 0.7-1.6" of the transit host stars, and are fainter by gen- 
erally 4 mag in the i band. Daemgen et al. found that their pres- 
ence changed the physical properties of the TEPs by roughly Icr. 
Of the three affected objects, TrES-2 and TrES-4 are analysed in 
the current work and WASP-2 is the subject of a separate publica- 
tion cSouthworth et al.201Q) . 



2.4.2 Accounting for third light 

The observations of Daemgen et al.' ('2009') make it possible to ac- 
count for third light when analysing transit light curves. They mea- 
sured magnitude differences (light ratios) in the SDSS i and 2 pass- 
bands. When necessary I have propagated the light ratios to other 
passbands by convolving synthetic spectra from ATLAS9 model at- 
mospheres l lKuruczill979l,ll993h with passband response functions 
made available by the Isaac Newton Grouflfl 

Armed with L3 values for the correct passbands, I have in- 
cluded these in the same way as e and uj. JKTEBOP was modified to 
accept measured L3 values as observations, and L3 was included as 
a fitted parameter. Note that several published studies have instead 
simply subtracted L3 from a light curve before analysis, which is 
stati stically incorrec t as it n eglects the uncertainty in L3 . 

IPaemgen et alj ( 12009 ) surveyed fourteen transiting systems 
and detected faint companions to three of them. The companions 



|http : / /cat server .ing.iac.es/filter/l 



3 INCORPORATING STELLAR MODELS : ABSDIM 

Analysis of a transit light curve gives the quantities ta, rb and 
From RV measurements of the parent star it is possible to obtain 
e, ui, and the velocity amplitude of the star, A'a- With these ob- 
servables we remain unfortunately one piece of information short 
of being able to calculate the full physical properties of the system. 
An additional constraint is needed, and this is generally supplied by 
forcing the properties of the star to match the predictions of theoret- 
ical stellar evolutionary models. To guide this process we can use 
the spectroscopically measured TcS and metal abundance, [^] , of 
the star. 

In the current work I adopt the method outlined in Paper II, in 
which the variable governing the solution process is taken to be Jfb, 
the orbital velocity amplitude of the planet. An initial value of Kh 



From transit light curves we also get Porb and Tq. The uncertainty in 
Porb is generally negligible, and To does not enter the ABSDIM analysis. 



© 0000 RAS, MNRAS 000, 000-000 



Homogeneous studies of transiting extrasolar planets. III. 



5 



Table 1. Physical ingredients and coverage of the stellar models used in this work. Yini is the primordial helium abundance, Ay/AZ is the helium-to-metals 
enrichment ratio, Zq is the solar metal abundance (fraction by mass) and Qmlt is the mixing length parameter. 



Model set 


Reference 


Range in 
mass ( M0) 


Range in metal 
abundance (Z) 


lini 


Ay 

AZ 




Q^MLT 


Notes 


Claret 


Claret (2004b. 2005. 2006. 2007) 


0.2 to 1.5 


0.01 to 0.05 


0.24 


2.0 


0.02 


1.68 


Calculated on request 


Y2 


Demarque et al. (2004) 


0.4 to 5.2 


0.00001 to 0.08 


0.23 


2.0 


0.02 


1.743 




Teramo 


Pietrinferni et al. (2004) 


0.5 to 10.0 


0.0001 to 0.04 


0.245 


~1.4 


0.0198 


1.913 




VRSS 


VandenBerg et al. (2006) 


0.4 to 4.0 


0.005 to 0.050 


0.23544 


2.2 


0.188 


1.90 




DSEP 


Dotter et al. (200|) 


0.1 to 5.0 


0.000041 to 0.0404 


0.245 


1.6 


0.0189 


1.938 





is defined, usually in the region of 150 km s ^ , and the full physi- 
cal p roperties of the system are calculated using standard formulae 
(e.g. lHilditch|[200lh . Armed with the resulting stellar mass, Ma, 
and [§] , I linearly interpolate within tabulated theoretical model 
results to find the predicted radius and T^s of the star. This process 
is iteratively repeated whilst varying Kh in order to minimise the 
figure of merit 



fom : 



(obs) 



+ 



(obs) 



(prcd) 



cff 



(3) 



which results in the best-fitting system properties. In principle it is 
possible to also solve for the age of the system, but in practise the 
wide variety of evolutionary timescales of stars make this difficult. 
I therefore step through the possible ages of the star in 0.1 Gyr in- 
crements, starting at zero age and finishing when the star leaves the 
main sequence, in order to find the best overall solution. I do not 
make any attempt to match spectroscopically measured log g val- 
ues as they are usually much less reliable than the surface gravity of 
the star calculated from the Ma and Ra obtained above. The above 
procedure implicitly applies the strong c onstraint on stellar density 
obtain ed from the light curve analysis dSeager & Mallen-OmelasI 
l2003t) . 

The uncertainties in the system properties are calculated by a 
perturbation analysis, in which each input parameter is modified by 
its la uncertainty and new solutions specified. The uncertainty for 
each output parameter is then calculated by adding the uncertainties 
due to each input parameter in quadrature. This perturbation analy- 
sis has the advantage of yielding detailed error budgets, where the 
effect of the uncertainty of every input parameter on every output 
parameter is specified. These error budgets indicate what additional 
observations are the best for improving our understanding of a spe- 
cific TEP. 



3.1 Which stellar models to use? 

As outlined above, the physical properties of TEPs are calculated 
by forcing the properties of the parent star to match theoretical ex- 
pectations. This dependence on theoretical predictions is a concern 
and will cause a systematic error. It is well known that whilst theo- 
retical models are pretty good at reproducing the actual properties 
of stars, the various model sets are not flawless and do not agree 
perfectly with each other. 

The existence of different theoretical model sets for low-mass 
stars opens the possibility of using several of them and explicitly 
deducing the systematic errors in TEP properties caused by their 
use. In Paper II six different sets of theoretical models were inves- 
tigated and three adopted for calculating the planet properties. The 
Siess and Cambridge-2007 models were in relatively poor agree- 



ment with other models, and the Cambridge-2000 models had a 
lower coverage of the relevant parameter space than other models. 
I was therefore left with only three different model sets, which was 
insufficient to define high-quality systematic error estimates. On 
top of this, the Padova models included heavy element abundances 
only up to Z = 0.03 so did not cover quite a few transiting systems. 

In the current paper I have therefore adopted the same solu- 
tion procedure as introduced in Paper II, but with a significantly re- 
vised database of theoretical model predictions and with one more 
change. Instead of using the Claret models to define my baseline 
solutions and two other model sets to obtain systematic error es- 
timates, I have used the unweighted mean and standard deviation 
of the results from all five model sets to describe the baseline solu- 
tions and systematic errors. The dependence of the final results on a 
single model set is therefore broken: all five model sets are treated 
equally and the choice of which sets of models to use becomes less 
important. 

Of the sets of theoretical models included in Paper II, only 
the models survive unchanged here (see TablefTJ. The Claret 
models have been supplemented by additional calculations for 
higher metal abundances of Z = 0.06 and 0.07. The third model 
set used here is Teramclfl ( Pietrinferni et alj|2004) . and I selected 
the ones with moderate convective core overshooting (for masses 
>1.1 M0), the standard mass loss law (77 = 0.4) and normal el- 
emental abundances (scaled-solar, i.e. no enhancement of the a- 
elements). For the fourth model set I acqui red the Victoria- Regina 
(VRSS) modelfl jVandenBerg et alj |20()^) with scaled-solar ele- 
mental abundances. In these models the convective core overshoot- 
ing parameter depends on mass and is empirically calibrated. The 
fifth and final model set (DSEP) comes from Dartmouth Stellar 
Evolution Databas^3 ™d again comprises the calculations with 
scaled-solar elemental abundances. I selected those models which 
follow the standard helium-to-metal enrichment law. The DSEP 
models include a contraction to the zero-age main sequence which 
can take tens of Myr. 

In Fig.|4]l compare the predictions in mass and radius of the 
five sets of stellar models, for an age of 0.5 Gyr and for the adopted 
solar chemical composition (which differs between models). The 
models have a fairly good agreement in the mass-radius plane, par- 
ticularly near 1.0 Mq as they are calibrated on the Sun, but a wider 
variety in the mass-Tcff plane. Also shown in Fig.|4]are the masses, 
radii and T^b values of the sample of detached eclipsing binary star 



8 Obtained from the BaSTI database on 17/11/2009: 
| http : //a lbione . oa-teramo . inaf . it/| 

^ Obtained on 18/1 1/2009 from http : / / www . cadc-ccda . hia-iha . 
nrc-cnrc.gc.ca/ cvo/ community/ Vict or iaReginaModels/. 
l'^ Obtained on 18/1 1/2009 from http : / / stellar . dartmouth . edu/ 
'^models/index . html. 
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Mass (Mq) Mass (Mg) 

Figure 4. Mass-radius (left) and mass-Tcfj (right) plots showing the predictions of the five sets of stellar models adopted in this work (solid lines). The 
predictions are for an age of 0.5 Gyr (to minimise the effects of evolution to and from the ZAMS) and for the solar chemical composition (which varies 
between models). The measured masses, radii and T^fi values of a sample of detached eclipsing binaries (see Paper II) are shown for compaiison, using blue 
errorbars. The Sun is indicated by the usual 0; note that the predictions of the models do not pass through the solar values on these plots as the Sun is much 
older than 0.5 Gyr. 



systems constructed in Paper II. It can be seen that the disagreement 
with the models is much larger than the errorbars. This conclusion 
holds for all chemical compositions for which the five sets of mod- 
els are available. Similarly, adopting an age either before or beyond 
the main sequence can provide an agreement for individual eclips- 
ing binaries, but not for all of them simultaneously. 

Fig.|4]illustrates what can be expected for the variation of sys- 
tematic errors with mass. The models agree with each other very 
well in some regions, so systematics will be minimised, and less 
well at lower and higher masses, when systematics will be larger. 
The agreement between models is clearly much better than with the 
properties of well-studied eclipsing binaries. This means that using 
the five model sets will lead to only a lower limit on the system- 
atic errors in the properties of TEPs. A probable upper limit to the 
systematic effects can be obtained by calculating solutions with an 
eclipsing binary mass-radius relation instead of a stellar model set; 
this is applied below and discussed further in Paper II. 



3.2 Calculating the physical properties of transiting planets 

Using the method and theoretical stellar models outlined above, the 
mass, radius, surface gravity and mean density of the star (Ma, -Ra, 
loggTA, Pa) and of the planet (A/b, -Rb, ffb, ph) can be calculated. 
For each TEP I have obtained results for each of the five stellar 
model sets and also using an empirical mass-radius relation de- 
fined by the eclipsing binaries. The final result for each parameter 
is the unweighted mean of the five stellar-model results. The statis- 
tical errorbar is taken to be the largest one from these five solutions 
and the systematic errorbar is taken to be the standard deviation of 
the parameter values from the five solutions. I include the final Kh 
value when possible to aid the comparison between different solu- 



tions for the same planet. K\, is the parameter through which all of 
the stellar model dependence enters. 

The surface gravity of the planet, can be calc ulated from 
purely geometrical observed quantities jSouthworth et al. 2007^ 
so has no systematic error. Similarly , the stellar density, pa, is in- 
dependent of the stellar models tSeager & Mallen-OmelaslFlOOBl) if 
Ma ^ Afb is assumed. 

In addition to the above parameters, I have also calculated the 
equilibrium temperature (Tcq) and Safronovi ( i l972i) number (0) of 
the planet. Tcq is given by the equation 

— (^)"'(^)'" 

where A is the Bond albedo and F is a heat redistribution factor. 
Because A and F are not known precisely I instead calculate a 
modified equilibrium temperature (T^'q) which equals Tcq \f A = 
1-4F: 




The Safronov number is defined as the ratio of the escape velocity 
to the orbital velocity of the planet: 

1 /Ks^y ^ /Mb\ ^Mb 

2VVo,b/ \R^,)\Ma) t^Ma 

From Eq.|5] above it can be seen that T^'q depends only on the 
stellar Tcfr and the fractional radius obtained from the light curve 
analysis. Tjq therefore turns out (like gb) to be independent of the 
stellar models, but does have some systematic error as it is depen- 
dent on the effective temperature scale of low-mass (F, G and K) 
stars. 
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Table 2. Parameters from the light curve analyses presented here and in previous works, and used here to determine the physical properties of the TEPs. The 
orbital periods are taken from the literature, and the bracketed numbers represent the uncertainty in the preceding digits. Systems for which orbital eccentricity 
was accounted for are indicated with a * in the column marked "e?". 



System 


Orbital period 


e? 


Orbital inclination. 


Fractional stellar 


Fractional planetary 


Reference 




(days) 




i (degrees) 


radius, r\ 


radius, 




GJ436 


2.64389524 (76) 




86.43 ±0.18 


0.0731 ±0.0027 


0.00605 ± 0.00023 


Paper I 


HAT-P-1 


4.4652934 (93) 




86.25 ±0.22 


0.0935 ± 0.0025 


0.01051 ±0.00031 


This work 


HAT-P-2 


5.6334729 (61) 


* 


85.9 ± 1.5 


0.1247 ±0.0106 


0.00847 ± 0.00082 


This work 


HD 149026 


2.8758911 (25) 




88.0 ±2.0 


„ ^ , 1-0 019 




Paper I 


HD 189733 


2.21857578 (80) 




85.78 ±0.25 


0.1113±0.0031 


0.0175 ±0.0005 


Paper I 


HD 209458 


3.52474859 (38) 




86.590 ± 0.046 


0.1 1384 ±0.00041 


0.01389 ±0.00006 


Paper I 


OGLE-TR-10 


3.101278 (4) 




83.87 ±0.69 


0.157 ±0.009 


0.0182 ±0.0011 


Paper I 


OGLE-TR-56 


1.211909 (1) 




79.8 ±2.4 


0.245 ± 0.026 


0.0241 ±0.0034 


Paper I 


OGLE-TR-111 


4.0144479 (41) 




88.11 ±0.66 


0.0842 ±0.0038 


0.01 107 ±0.00067 


Paper I 


OGLE- IR- 113 


1.4324757 (13) 




87.7 ±1.4 


0.1592 ± 0.0043 


0.02331 ± 0.00089 


This work 


OGLE-TR-132 


1.689868 (3) 




83.3 ±2.4 


0.211 ±0.020 


0.0198 ±0.0024 


Paper I 


OGLE-TR-182 


3.97910(1) 




84.3 ±1.2 


0.137 ±0.014 


0.0135 ±0.0013 


This work 


OGLE-lR-21 1 


3.67724 (3) 




88.0 ± 2.0 


n 1/199+0.0150 


n n-T 1 01 +0.00146 
U.UiiSl_o 00083 


This work 


OGLE-TR-L9 


2.485533 (7) 




82.07 ± 0.69 


0.1731 ±0.0083 


0.01910 ±0.00085 


This work 


TrES-1 


3.0300728 (6) 




88.67 ±0.71 


0.0964 ±0.0018 


0.01331 ±0.00035 


Paper I 


TrES-2 


2.4706101 (18) 




83.80 ±0.36 


0.1282 ±0.0035 


0.01658 ± 0.00043 


This work 


TrES-3 


1.3061864(5) 




82.07 ±0.17 


1 n noi 7 

n ifififi+ 


n n97'^ 1 """^"^^^^^ 

U.UZ/0±_Q QQQ43 


This work 


TrES-4 


3.553945 (/5) 




o 1 CT 1 n 


n 1 om 1 n aaoi 

0.1802 ± 0.0083 


r\ A 1 "7 /I 1 A AA 1 n 

U.Ol /4 ± U.U012 


Ihis work 


WASP-1 


2.519961 (18) 




88.0 ±2.0 


n 1 7Q7+U.UU57 


0.0182lO;°«?^ 


Paper I 


WASP-2 


2.15222144 (39) 




84.81 ±0.17 


0.1238 ±0.0018 


AA1/^^T 1 AAAAOA 

0.01643 ± 0.00030 


C _4-l_ J- 1 /'^Al A\ 

Southworth et al. (2010) 


WASP-3 


1.846835 (2) 




84.1 ± 1.3 


0.201 ±0.010 


0.0218 ±0.0011 


This work 


WASP-4 


1.33823150 (61) 




89.0 ± 1.0 




n 09819+0-00022 
U-U^812_o 00014 


Southworth et al. f2009b) 


WASP-5 


1.6284246 (13) 




85.8 ± 1.1 


0.1847 ±0.0061 


0.02050 ±0.00091 


Southworth et al. f2009a) 


WASP- 10 


3.0927616 (10) 


★ 


88.81 ±0.40 


0.0865 ±0.0041 


0.01349 ±0.00065 


This work 


WASP- 18 


0.94145181 (44) 




85.0 ±2.1 


0.2795 ±0.0084 


0.0272 ±0.0012 


Southworth et al. (2009c) 


XO-1 


3.9415128 (28) 




89.06 ±0.84 


0.0886 ±0.0019 


0.01 166 ±0.00035 


Paper I 


XO-2 


2.6158640(16) 




88.8 ± 1.2 


n -,907+0.0024 

U--l-^o< _Q 0047 


U.UiOUU_Q 00070 


This work 


XO-3 


3.1915289 (32) 


* 


83.89 ±0.40 


0.1447 ±0.0046 


0.01317 ±0.00047 


This work 


XO-4 


4.12502 (2) 






n 1 Qnn+O-0283 
U.iciUU_Q 0051 


01124+0-00334 
U.Uilz^_0 00054 


This work 


XO-5 


4.187757 (11) 




87.04 ±0.65 


0.1004 ± 0.0049 


0.01054 ± 0.00073 


This work 



4 RESULTS FOR INDIVIDUAL SYSTEMS 

In this section I present the photometric (JKTEBOP) and absolute- 
dimensions (ABSDIM) analyses of a set of thirty TEPs based on 
high-quality data. In many cases I adopt the JK TE BOP re sults from 
Paper I or later works jSouthworth et al.l2OO9allbll3.l2O10l) . The final 
JKTEBOP results of all TEPs are collected in Tabled which also 
includes the orbital periods and indicates for which systems a non- 
circular orbit was adopted. The mass ratio of each TEP is required 
as an input parameter for the light curve analysis, but in all cases 
its effect on the solution is negligible. Representative values have 
been taken from the literature but will not be discussed further. 

The physical properties of all thirty TEPs are obtained or re- 
vised in the current work, using the new theoretical model sets 
discussed in Sect.[3] This also requires Toff, [^1 and Ka values 
for each system. These are summarised in Tal3le[3] The values are 
mostly unchanged for the fourteen TEPs studied in Paper II, but in 
a few cases improved values have become available and replace 
the previous entries. In Papers I and II the individual systems were 
tackled roughly in order of increasing complexity. The current work 
reverts to the more structured approach of attacking the TEPs in 
alphabetical order, beginning with those objects for which a light 
curve analysis is presented (Sects. [43] to 14. 15l then moving on to 
those whose photometric parameters are adopted unchanged from 
PaperI(Sect. l4J6t . 



4.1 HAT-P-1 

HAT-P-1 was found to be aTEP by'Sakos et al' ('2007a') from data 
taken by the HAT survey dSako s et al. 2002, 2004). Its low mass 
(0.5Mjup) and large radius (1.2Rjup) make it one of the least 
dense exoplanets known. Excellent light curves from the FLWO 

1.2 m (2 band), Lick 1.0 m Nickel (Z band) and Wise 1.0 m tele- 
scopes were presented bvl Winnetal] ( 12007 c') and the first two of 
these were included in Paper I. Since then additional data from the 
Nickel (Z band) and the Ha waiian 2 m Magnum {V band) tele- 
scopes have been obtained bv lJohnson et al.l i l2008l) . In this work I 
have analysed the latter two datasets in order to refine the results 
from Paper I. In both cases I have adopted solutions with the lin- 
ear LD coefficient fitted and the nonlinear coefficient fixed but per- 
turbed in the Monte Carlo simulations ('LD fit/fix'). The residual 
permutation analyses indicate that correlated errors are important 
for both datasets. 

The final light curve parameters are the weighted means of 
those for the four studied datasets. The results agree well with each 
other except for k, for which xS = 2.8. The errorbar for k has 
been multiplied by \/2.8 to account for this. The light curve fits 
are plotted in Fig.|5]and summarised in Table A3. They are in good 
agreement with literature values. 

The physical properties of HAT-P-1 have been calculated us- 
ing the five different sets of stellar evolutionary models plus the 
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Table 3. Measured quantities for the parent stars which were adopted in the analysis presented in this work. 



System 


Velocity amplitude ( m s ^ ) 


%B (K) 


Reference 






. H . 


Reference 


GJ436 


18.34 ±0.52 


Maness et al. (2007) 


3500 ± 100 


Bean et al. (2006) 


-0 


03 ± 





2 


Bonfils et al. (2005) 


HAT-P-1 


59.3 ± 1.4 


Johnson et al. (2008) 


5975 ± 50 


Bakos et al. (2007a) 





13± 





05 


Bakos et al. (2007a) 


HAT-P-2 


983.9 ± 17.2 


Pal et al. (2010) 


6290 ± 60 


Pal et al. (2010) 





14± 





08 


Pal et al. (2010) 


HD 149026 


43.3 ± 1.2 


Sato et al. (2005) 


6147 ±50 


Sato et al. (2005) 





36 ± 





05 


Sato et al. (2005) 


HD 189733 


200.56 ±0.88 


Boisse et al. (2009) 


5050 ± 50 


Bouchv et al. (2005b) 


-0 


03 ± 





05 


Bouchv et al. (2005b) 


HD 209458 

OGLE-TR-10 

OGLE-TR-56 


85.1 ± 1.0 
80 ± 17 
212 ±22 


Naef et al. (2004) 
Konacki et al. (2005) 
Bouchy et al. (2005a) 


6117±50 
6075 ± 86 
6119 ±62 


Santos et al. (2004) 
Santos et al. (2006) 
Santos et al. (2006) 







02 ± 
28 ± 
25 ± 







05 
10 
08 


Santos et al. (2004) 
Santos et al. (2006) 
Santos et al. (2006) 


OGLE-TR-111 


78± 14 


Pont et al. (2004) 


5044 ±83 


Santos et al. (2006) 





19± 





07 


Santos et al. (2006) 


OGLE-TR-113 


267 ± 34 


TWH08 


4804 ± 106 


Santos et al. (2006) 





15± 





10 


Santos et al. (2006) 


OGLE-TR-132 


167 ± 18 


Moutou et al. (2004) 


6210 ±59 


Gillon et al. (2007) 





37 ± 





07 


Gillon et al. (2007) 


OGLE-TR-182 


120 ± 17 


Pont et al. (2008) 


5924 ± 64 


Pont et al. (2008) 





37 ± 





08 


Pont et al. (2008) 


OGLE-TR-211 


82± 16 


Udalski et al. (2008) 


6325 ±91 


Udalski et al. (2008) 





11± 





10 


Udalski et al. (2008) 


OGLE-TR-L9 


510± 170 


Snellen et al. (2009) 


6933 ± 58 


Snellen et al. (2009) 


-0 


05 ± 





20 


Snellen et al. (2009) 


TrES-1 


115.2 ±6.2 


Alonso et al. (2004) 


5226 ± 50 


Santos et al. (200&) 





06 ± 





05 


Santos et al. (2006) 


TrES-2 


181.3±2.6 


O'Donovan et al. (2006) 


5795 ± 73 


Ammler-von Eiff et al. (2009) 





06 ± 





08 


Ammler-von Eiff et al. (2009) 


TrES-3 


369±11 


Sozzetti et al. (2009) 


5650 ±75 


Sozzetti et al. (200£) 


-0 


19 ± 





08 


Sozzetti et al. (2009) 


TrES-4 


97.4 ±7.2 


Mandushev et al. (2007) 


6200 ±75 


Sozzetti et al. (200^ 





14 ± 





09 


Sozzetti et al. (2009) 


WASP-1 


111±9 


Wheatley et al. (2010) 


61 10 ±50 


Stempels et al. (2007) 





23 ± 





08 


Stempels et al. (2007) 


WASP-2 


153.6 ±3.0 


Triaud et al. (2010) 


5 150 ±80 


Triaud et al. (2010) 


-0 


08 ± 





08 


Triaud et al. (2010) 


WASP-3 


zyu. J ± y. J 


iripatnl et al. (zUlU) 


6400 ± 100 


Pollacco et al. (2008) 





00 ± 


U 




roiiacco et al. (ZUUJS) 


WASP-4 


242.113^ 


Triaud et al. (2010) 


5500 ± 100 


Gillon et al. (2009) 


-0 


03 ± 





09 


Gillon et al. (2009) 


WASP-5 


268.7 ± 1.8 


Triaud et al. (2010) 


J /UU ± lUU 


/"'lu^.T r^t nl /oonn\ 
Gillon et at. (zUUy) 





09 ± 





09 


Gillon et al. (2009) 


WASP- 10 


553.1 ±7.5 


Johnson et al. (2009b) 


4675 ± 100 


Christian et al. (2009) 





03 ± 





20 


Christian et al. (2009) 


WASP- 18 


1816.9 ±2.0 


Triaud et al. (2010) 


6400 ± 100 


Hellier et al. (2009) 





00 ± 





09 


Hellier et al. (2009) 


XO-1 


116.0 ±9.0 


McCuUough et al. (2006) 


5750 ±50 


McCulloush et al. (2006) 





02 ± 





05 


McCullough et al. (2006) 


XO-2 


85 ±8 


Burke et al. (2007) 


5340 ± 50 


Burke et al. (2007) 





45 ± 





05 


Jurke et al. (2007) 


XO-3 


1488 ± 10 


Winn et al. (2009^ 


6429 ± 75 


Johns-Krull et al. (2008) 


-0 


1 


8± 





05 


'ohns-Krull et al. (2008) 


XO-4 


163 ± 16 


McCulloush et al. (2008) 


6397 ± 70 


McCulloush et al. (2008) 


-0 


04 ± 





05 


McCullough et al. (2008) 


XO-5 


144.9 ±2.0 


Pal et al. (2009) " 


5370 ±70 


Pal et al. (2009) 





05 ± 





06 


Pal et al. (2009) 



empirical mass-radius relation from Paper II. The individual solu- 
tions are given in Table A4 and then compared with literature val- 
ues, where a good agreement is found. 



4.2 HAT-P-2 



HAT-P- 2 was discovered to be a TEP system by I Sakos et"aD 
( l2007bl) . under the name HD 147506. It is a very bright system 
(V = 8.7) with a massive planet (Mb = 8.74 Mjup) in a highly 
eccentric o rbit. It has been found not to exhibit a spin-orbit mis- 
alignment jWinn et al.ll2007j : iLoeillet et al.ll2008h. in contr ast to 
other massive TEPs on eccentric orbits ('J ohnson etalj200 9aV 

Good 2:-band light curves of HAT-P-2 have been published by 

iBakos eta l. (2007b), covering one transit with the FLWO 1.2m, 



and by iPal et al.. t201Q) , covering another six transits with the 
FLWO 1.2 m and the Perkins telescopes. Here I analyse the FLWO 
datasets together, omitting the small amount of data taken on 
the night of 2007/03/18, as well as the Perkins data. One com- 
plication is the orbital eccentricity: this was accounted for us- 
ing the method discussed in Sect. |2.3| and adopting the constraints 
ecoso; = -0. 5152 ± 0.0036 and esinw = -0.0441 ± 0.0084 
jPal et al.l2oiob . In both cases correlated errors were not important 
and the LD fit/fix solutions were adopted. The best fits are shown 
in Fig.[6] 

The two light curve solutions unfortunately do not agree very 
well (9.3cr discrepancy in k). I therefore adopt the FLWO 1.2 m 
results, as these are the much more extensive of the two sets of 
data and have full coverage of the transit phases. The FLWO results 



agree well with those o f lPaletal-lJioiB ), for which most of the data 
come from, but have a larger ta and rb than other literature values 
(Table A7). 

As expected given the light curve results, my ABSDI M analysis 
return s system properties in good agreement with those of |Pal et al] 
1L2OIO) but not with other literature studies (Table A8). The prime 
mover in the most recent solutions is ta, which has a strong effect 
on the density of the star and thus the other physical properties. The 
radius of the planet is uncertain by 10%, despite the existence of a 
high-quality light curve for HAT-P-2, because the transit depth is 
shallow (0.6%). An improved photometric study is warranted. 



4.3 OGLE-TR-113 

Like OGLE-TR-132 (studied in Paper I), OGLE-TR-1 13 w as iden- 

'lUdalski et al] ( l2002h and its 



tified as a possible plane tary system by 



nature was confirmed bv lBouchv erai] l l2004 using the OGLE light 



curve and n ew RV measurements . It was independently confirmed 
as a TEP bv lKonacki et alj ( l2004i) . also from the OGLE light curve 
and high-precision RVs, a nd an abundance a nalysis of the parent 
star has been presented bv lSantos et~ai] ( l2006h . Whilst OGLE-TR- 
113 exhibits a deep transit, its photometric tractability is hindered 

by the presence of a brighter star only 3" away. 

Apart from the OGLE discovery observations jUdalski et al.l 

l2002h . three photometric s tudies of OGLE-TR-113 have been pub- 
lished. iGillon et alj ( I2OO6I) used the ESO New Technology Tele- 
scope (NTT) and SUSI2 imager to observe two transits in the R 
band, and obtained what is currently the best light curve of OGLE- 
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Figure 5. Phased light curves of HAT-P- 1 compared to the best fits found 
using JKTEBOP and the quadratic LD law in Paper I and in this work. The 
best fits and residuals are offset from unity and zero fluxes, respectively, 
for di splay purposes. Th e light curves a re, from top to bott om, FLWO z- 
band jwinn et al.ll2007cl) Lick Z-band iwinn et alj|2007cl). Magn um V- 
band jjohnson et al.l2008l) and Nickel Z-band jjohnson et alj2008l) . 



1.00 




0.98 



0.96 



0.94 




-0.02 



-0.01 0.00 0.01 
Orbital phase 



0.02 



Figure 6. Phased z-band light curves of HAT-P-2 compared to the best fits 
found using JKTEBOP and the quadratic LD law. The best fits and residuals 
are of fset for display purposes. The upp er light curve is from the FLWO 



1.2m 
1.8m 



set tor display purposes. Ihe upp e 
Bakos et ia.l2007b| - |Pal et alj201Cll) 



Pal et al j2O10|y 



and the lower is from the Perkins 



a few Gyr, and in several cases is up against the edge of the stellar 
model grid at 20 Gyr Aside from that, the properties of the star 
and planet are rather well-determined but would benefit from an 
improved Ka value as well as a better light curve. The agreement 
with literature studies is good, although it seems that in some cases 
the published errorbars are smaller than one would expect. 



TR-1 13. [Snellen & Covinol ilOOH) observed a JsT-band transit and 
an occultation of the system using NT T/SOFI, and detect ed the oc- 
cultation with a significance of 2.Ha. IPfaz et al] ( l2007t) obtained 
l/-band photometry of one transit using a Very Large Telescope 
(VLT) and the VIMOS instrument; additional data taken in the I 
and Ks bands are unavailable and of lower qual ity. 

In this work I analyse the iGillon et al.l observations, the 

[Snelle n & Covino transit light curve, and the l/-band data obtained 
by IPiaz et alj For the second of these three datasets I allowed for 
light from the planet with a surface brightness ratio of 0.07 ± 0.02. 
The surface brightness ratio is a parameter of the JKTEBOP model 
which is important for eclipsing binary systems but usually left at 
zero for transiting systems due to the faintness of the planet with 
respect to the star. The best fits are shown in Fig.|7]and given in 
Table A12. In all three cases correlated noise is not important. For 
the Snellen light curve I had to adopt solutions with both LD coeffi- 
cients fixed, but for the other two I was able to use the LD fit/fix so- 
lutions. The final results for the Gillon and Snellen data are in good 
agreement. The solution of the Diaz data prefers a rather higher i 
and lower ta and Tb. I therefore combined the Gillon and Snellen 
data results to obtain the final light curve parameters. 

The resulting physical properties of OGLE-TR-113 are given 
in Table A 13. The system age is constrained only to be more than 



4.4 OGLE-TR-182 

OGLE-TR-182 is the sixth TEP discovered as a result of the OGLE 
search for light variability in selected fields in the Sout hern hemi- 
sphere. Its discovery and analysis was presented by IPont et al] 
OOOS) , which remains the only study of this object to date. OGLE- 
TR-182 is difficult becau se of its faintness (V — 16.8 and / — 
15.9) and crowded field. ' Pont et al.l feOOSi) obtained 24 RV mea- 
surements using VLT/FLAMES/UVES, and a light curve with 
VLT/FORSl. 

The VLT light curve is analysed here and is rather affected 
by correlated noise. Including the linear LD coefficients as a fitted 
parameter gives substantially better fits than with both LD coeffi- 
cients fixed, but the data cannot support the determination of both 
LD coefficients. I therefore adopt the LD fit/fix solutions (seeFigJSl 
and Table A 15), which are not in good agreement with lPont et al.l 
(2008). Compared to these authors I find a solution with a lower i 
and a correspondingly larger star and planet. 

The physical properties of OGLE-TR-182 are summarised 
in Table Al 6 and point to a planet with a rather low density of 
0.33 pjup. Howe ver, my results are rather different to those of 
iPont et al. liooi), and are in poorer agreement with the measured 
spectroscopic Tofi and (rather uncertain) log g measurement. My 
analysis procedure is more sensitive to the quality of the light curve 
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Figure 9. Phased VLT light curve of OGLE-TR-211 from ludalski et alj 
( 2008) compared to the best fit found using JKTEBOP and the quadratic LD 
law. The residuals are offset from zero for display purposes. 
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Figure 7. Phased light curves of OGLE-TR- 1 1 3 with the best fits found us- 
ing JKTEBOP and the quadra tic LD law and residuals of t he fits. From top 
to bottom the datasets aie Gillon et al U2OO6) (i?-band), Snellen & Coving 
2007) (Xs-band, binned by x5 for display purposes) and iDfazeTaT 



200711 (V-band). 



.00 



2 0.99 




0.98 -. 



0.97 - 



.-i ' Ji^/- ' :.%V 



-0.04 



-0.02 



0.00 
Orbital phase 



0.02 



0.04 



Figure 8. Phased VLT light curve of OGLE-TR- 1 82 from lPont etaP i2008h 
compared to the best fit found using JKTEBOP and the quadratic LD law. 
The residuals are offset from zero for display purposes. 



than the more 'global' approach taken by Pont et al., so is poten- 
tially more susceptible to correlated noise. This possibility should 
be investigated by acquiring a new light curve, under good seeing 
conditions to cope with the crowded field. 



4.5 OGLE-TR-211 

OGLE-TR-211 is th e seventh TEP discovered using OGLE data 
dUdalski et ai]|2008|) . Its relative faintness means that the available 
follow-up photometry and spectroscopy is not definitive. The par- 
ent star is more massive and also more evolved than the Sun, which 
results in the transit being rather s hallow (Fig.[9] >. Here 1 analyse 
the VLT light curve presented by Udalski et al. ( 2008), ignoring 
the observational errors supplied with the data which are quite dis- 
cretised (the only values are 0.001, 0.002 and 0.003) and contribute 
to instability in the light curve solution. 

1 am not able to get a determinate solution to the VLT data. 
Possible fits occupy a locus extending from a high i with low ta to 
a lower i with a large ta- I have therefore calculated solutions for 
a range of i values and retained only those which in the ABSDIM 
analysis result in a Tcft within a conservative 3cr of the observed 
value. The observed stellar log g did not provide a useful constraint. 
Allowable solutions extend from i = 90° down to a sharp cut-off 
around i — 86.25° so I present solutions for i = 86, 88, and 90° 
in Table A17. For the final result I accept the LD fit/fix solutions 
for i — 88° but specify errors which account for both the Monte 
Carlo errorbars and the variation between the different solutions 
(Table A 18). The correlated errors are again important, as can be 
seen in Fig.|9] 

The physical properties of OGLE-TR-211 are shown in 
Table A19 and are in reasonable agreement with those of 
lUdalski et alj ( l2008h except for the planetary mass. A/b depends 
mainly on the measured Ka, for which both studies have used the 
same value, so it is not clear why such a discrepancy should arise. 
Table A19 includes the first determinations of the age and density 
of the star, planetary equilibrium temperature (which is quite high 
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Figu re 11. Phased light curve of the transits of TrES-2 from lHolman et all 
(2007) compared to the best fit found using JKTEBOP and the quadratic 
LD law with the linear LD coefficient included as a fitted parameter The 
residuals of the fit are plotted at the bottom of the figure, offset from zero. 
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Figure 10. Phase d GROND light curves of OGLE-TR-L9 from 
ISnellen et all i2009l) compared to the best fit found using JKTEBOP and the 
quadratic LD law. From top to bottom the light curves are g, r, i and z. The 
residuals are offset from zero to the base of the figure. 



relatively well determined (2.6^0 7 Iq * Gyr) because the star has 
evolved away from the zero-age main sequence. OGLE-TR-211 
would certainly benefit from additional spectroscopic and photo- 
metric observations. 



4.6 OGLE-TR-L9 

OGLE-TR-L9 was disc overed within the O GLE-II survey data 
( IUdalskietai]|l997h by ISnellen et al] j2009l) . and is a relatively 
massive planet orbiting a rapidly-rotating (Vsini = 39kms^^) 
F3V star. High-quality follow-up light curves were obtained by 
Snellen et al. using the newly-commissioned GROND instrument 
( iGreiner et all [2008) on the 2.2 m telescope at ESO La Silla. 
GROND is a CCD imager which utilises dichroics to observe si- 
multaneously in seven passbands (SDSS griz and near-infrared 
JHK). In the case of OGLE-TR-L9 the JHK data were too noisy 
to be useful, but the griz data are of good quality (Fig.llOl). 

The griz observations have been analysed here (Tables A20 to 
A23). The gri light curves are good enough to support LD fit/fix so- 
lutions but for the z data both LD coefficients were held fixed. The 
parameters for the four light curves were combined to obtain the 
final pho tometri c res ults (Table A24). Their agreement with those 
of Snellen et al. I I2OO9,) is not good - i and ta are correlated and 



my solution corresponds to a significantly higher i (3a) and lower 
rA and (both 6a). 

The ABSDIM results for GLE-TR-L9 (Table A 25) are in rea- 
sonable agreement with those o f lSnellenetaLl ( l2009l) , which is sur- 
prising given the differences in the photometric parameters. One 
reason for this is that the ABSDIM solution is governed mainly by 
the observed TcH (which has a relative uncertainty of 0.8%) rather 
than by ta (5%). I find that 0GLE-TR-L9b has one of the high- 
est rJqS (2039 ± 51 K) of any known TER Aside from its faint- 
ness, this planet is an excellent candidate for multicolour 'transmis- 
sion photometry' to detect variations in planetary radius with w ave- 
length due to atmospheric opacity effects ( Fortnev et al .I2OO8I) . The 
system would benefit from additional spectroscopy to provide im- 
proved measurements of [^] and Ka- This subsection completes 
my analysis of the TEPs discovered from OGLE survey data. 



4.7 TrES-2 

TrES-2 was th e second TEP discovere d by the Trans-Atlantic Exo- 
planet Survey jO'Donovan et al.ll2006l) and is both larger and more 
massive than Jupiter. Excellent grou nd-based light curve s of three 
transits were obtai ned and studied by Holman et ai]l l2007h , in close 
collaboration with lSozzetti et al' I( l2007h . The relatively low orbital 
inclination of this system (i = 83.8°) help ta and r^ to be deter- 
mined to a high accuracy. 

The Holman z-band light c urve of TrES-2 was stu died in Pa- 
per L but is revisited here because lOaemgen et al J ( 12009 !) have since 
found a fainter star separated from TrES-2 by 1.089 ± 0.008 arc- 
sec and with a magnitude difference of Az = 3.429 ± 0.010 (see 
Sect.|231l. 

TrES-2 is of additional interest because i Mislis & SchmittI 
(2009() found a possible decrease in the transit duration between 
their own and Holman's observations. This would most likely in- 
dicate that the orbital inclination is getting lower, which in turn 
point s to the pres ence of a low-mass third body in the TrES-2 sys- 
tem. [Scuderietalj (.200^) obtained new data which did not confirm 
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this hypothesis, but lMishs et al. I (I2OO9I) have since resurrected the 
changing i. In addition. iRabus et aT presented a transit tim- 

ing study of TrES-2 which found a small sinusoidal perturbation 
with a 0.2 day period, but with only moderate statistical signifi- 
cance. TrES-2 is i n the field of view of the NASA Kepler satellite 
dKoch et alj|2010h . so a light curve of remarkable quality should 

become available for future investigations of these possibilities^ 

Here I reanalyse the z-band light curve of Holman et al.l 

( l2007h . this time with the incorporation of a third light value of 
L3{z) = 0.0408 ± 0.0004. From the adopted T ^s of TrES-2 A 
(Tabl ed and the magnitude differences in i and z l lDaemgen et al.l 
I2OO9I) I find that the fainter companion star has Tes = 4390 ± 70 K 
(cons istent with the ~K5 spectral type found by lOaemgen et al.] 
l2009h . and is substantially further away than TrES-2 so is not phys- 
ically bound to the system. 

The results of the JKTEBOP analysis are given in Table A26 
and I adopt the LD fit/fix results. Correlated noise is not impor- 
tant. Table A27 shows a comparison with the results from Paper I, 
in which the analysis did not account for the faint companion star: 
k has decreased by la whilst ta + r^, and ta become smaller by 
less than 0.5ct. The best fit is shown in Fig.ll II 

The physical properties of the TrES-2 system are summarised 
in Table A28. My results agree with literature studies within the er- 
rors, although AIa and A/b are larger by roughly la. This can be 
traced back to the slightly smaller ta found above. More precise 
spectroscopic Tcs and [^] values would allow improved system 
parameters to be obtained, as would better photometry. This is cur- 
rently being obtained by the Kepler satellit e, and a stunning ligh t 
curve of TrES-2 can already be inspected in lGilliland et alj bOlOt) . 

4.8 TrES-3 

TrES-3 is one of the more massive TEPs (Mb = 1.91 Mjup) and 
orbits a rather cool star (Tch = 5650 K) with an orbital period of 
only 1.3 d. It was identified as a TEP by lO'Donovan et alj j2007h 
and a discovery-qua lity light curve has also been obtained by the 
SuperWASP survey JCoUier Cameron et alj2007[). FoUow-up tran- 
sit photomet r y has been presented by Sozze tti et alj ( l2009h and 
iGibson et alj |2009|) . and occultation (second ary eclipse) observa- 
tions have been secured bv lWinn et alj i2008a[). A space-bas ed light 
curve was obtained by the EPOXI mission iBallard et alj2009f) but 
has not y et been published. 

The lO'Donovan et al. B- and 2:-band photometry of 

TrES-3 is o f good quality, and the transit observation s of the follow- 
up papers dSozzetti et ^|2009| ; [oibson et alj|2009h are excellent. 
This, plus the fact that TrES-3 has a relatively low orbital inclina- 
tion (82°), means that the light curve parameters can be obtained 
to an unusually hig h precision. In this wo rk I analyse the B- and 
a-band data f rom lO'Dono van et alj ( l2007h . the V gri-hw\d obser- 
vations from I Sozzetti et alj 1 2009h. and the Liver pool Telescope 
(LT) RISE measurements from lGibson et al.l | |2009|) . The last of the 
datasets, which has a passband of approximately V+R, was sorted 
in phase and then binned down by a factor of 20 (from 1 1 350 to 
568 datapoints) to ease the computational burden. 

For the solutions of the seven light curves I adopt those with 
both LD coefficients fixed for V , r and i, and the LD fit/fix solu- 
tions for the remainder (Fig.ll2t. Correlated noise is important for 
the V , g and LT datasets. For most of the datasets I find that the 
residual-permutation algorithm returns errorbars which are signifi- 
cantly asymmetric, with larger upper errorbars than lower errorbars 
for both ta and rb, but that the Monte Carlo errorbars are close to 
symmetric. This implies that the red noise in the light curves is 
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Figure 12. Phased light curves of the transits of TrES-3 compared to 
the best fits found using JKTEBOP and the quadratic LD law. Successive 
datasets and residuals are offset in flux for display puiposes. From top to 
bottom are the B an d z data (O'Donovan et al. 2007 ), Vqri observat ions 
JSozzetti et al.l2009l) and LT V+i? measurements iGibsonet al .1120091) . 



biasing the results towards larger component radii, and would not 
have come to light if correlated noise was accounted for simply by 
rescaling the observational errors. 

I accordingly end up with asymmetric errorbars for each light 
curve solution. The V and g solutions are quite uncertain and dis- 
crepant with the other results, so were rejected. To calculate the 
final photometric result I combined the solutions of the remaining 
five light curves by multiplying their probability density functions. 
The result is given in Table A36 and is found to be in good agree- 
ment with literature results. 

The physical properties of the TrES-3 system are summarised 
in Table A37, and agree well with published values. For all five 
stellar model grids the best solution was found for zero age, and 
there is a possibility that edge effects will cause the uncertainties 
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Figure 13. Phased light curves of the transits of TrES-4 compared to the 
best fits found using JKTEBOP and the quadratic LD law. The upper data are 
in the z band and the lower data are in the B band. The fits and residuals 
are offset in flux for display purposes. 
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to be slightly underestimated for this object (for an example of this 
phenomenon see HD 189733 in Paper I). 



4.9 TrES-4 

The p lanetary system TrES-4 was discovered bv lMandushev et aP 
and is noteworthy for having a very hot (Tj^ = 1861 K) 
and low-density (0.15 pjup) planet. Revised physi cal properties 
of the system have been presented by TWH08 and lSozzetti et alj 
( |2009|) . and near-IR observations of the secondary ec lipse show no 
evidence for orbital eccentricity jKnutson et alj2009h . 

Lik e TrES-2 and W ASP-2, high-resolution imaging observa- 
tions bv iDaemgen et alj | [2009.) have detected a faint companion to 
TrES-4. It resides at an angular separation of 1.555 ± 0.005 arcsec 
and is fainter by Ai = 4.560 ± 0.017 and Az = 4.232 ± 0.025 
than TrES-4 A. These result in third light contributions of ^3(1) — 
0.0150 ± 0.0002 and L3{z) = 0.0199 ± 0.0005, which are taken 
into account in the light curve analysis as described in Sect. |2.4| 

In this work I study the hig h-precision B and z lig ht curves 
presented in the discovery paper (^ Mandushev et al.l2007l) . and pro- 
vide the first results to fully account for the third light contribu- 
tion. Using ATLAS9 model spectra I have propagated the third 
light to the B-band, finding LaiB) = 0.0040 ± 0.0003 and a 
Teff = 4206 ± 78 K for the companion star. The fainter star is 
more than twice as distant as TrES-4, so is not bound to the plane- 
tary system. 

The B data are rather sparse and do not allow LD coefficients 
to be fitted for (Fig.ll3t. The z observations cover two transits so 
were studied with Porb left as a fitted parameter; the LD fit/fix 
solutions were adopted. Correlated errors are important for both 
datasets. The two light curve solutions agree reasonably well (1.5cr) 



Figure 14. Phased light curves of WASP-3 compai'ed to the best fits found 
using JKTEBOP and t he quadratic LD law. From top to bottom the light 
curves are LT V+R taibsonet all [2008). whi ch have been binned by a 
factor of ten for analysis, Keele R and lAC / jPollacco et ai]|2008l) . The 
residuals are offset from zero to the base of the figure. 

within the rather large errorbars, and were combined by weighted 
mean to find the final parameter values. This solution corresponds 
to a smaller i, and larger star and planet than found previously (Ta- 
ble A40). From Fig.[T]we would expect that accounting for third 
light would lead to a smaller ta and larger rb and i, which is only 
partially concordant with the current situation. The variation in pa- 
rameter values must therefore be due to the different analysis meth- 
ods used. 

The physical properties of TrES-4 are given in Table A41 and 
conform to a more evolved star than found by previous studies, as 
expected for the larger ta found above. The properties of the planet 
agree well with literature values, but are quite uncertain. The AB- 
SDIM error budget shows that an improved light curve and addi- 
tional RV measurements would benefit this system. All four TEPs 
discovered by TrES have now been analysed in the current series 
of papers. 

4.10 WASP-3 

WASP-3 was identified as a possible TEP by 'Street et alj ( l2007h 
based on SuperWASP data ( PoUacco ^t al. 2006). Co nfirma tion of 
its planetary nature was provided by [Pollacco et al] ( l2008h . who 
presented four transit light curves obtained from various sources 
and found WASP-3 b to be one of the most strongly irradiated TEPs 
(Tcq — 2028 K). High-precision lig ht curves from the LT have 
since been presented and analysed bv loibson et alj ( |2008[) . Here I 
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study these data , plus the Keele fl-ban d and lAU 80 cm /-band ob- 
servations from IPoUacco et ah The other two hght curves 
(lAC 80 cm V and SuperWASP) have either large systematics or a 
large scatter. I have binned consecutive sets of 10 datapoints of the 
LT light curve in order to limit CPU time; the sampling rate of the 
binned data is 30 s. 

For the LT data I adopt the LD fit/fix results and find signifi- 
cantly larger errorbars from the residual-permutation analysis; cor- 
related noise is clearly visible in these data in Fig.[T4] The other 
two datasets also contain significant red noise, and the LD fixed so- 
lutions were adopted. The agreement between the three light curve 
solutions is excellent so they have been combined into a weighted 
mean (Table A45). Apart from a 2a large r rb, th e fina l values agree 



well with the studies of IPollacco et alJ ( l2008h and Gibson et al 



2009) 



2008 ) but not with the preliminary results given bv lDamasso et al 



The physical properties of the WASP-3 s ystem were originall y 
calculated using Ka = 251.2 ± 9.3 km s"^ (pollacco et ai]|2008h . 
After this w ork was comple t ed a re vised Ka of 276.0 ± 11.0 was 
presented by Simpson et alj | |2009|) . who also found that the angle 
between the planetary orbit and the stellar spin was A — 15lg" °. 
Shortly before the curre nt work was submitt ed a further study of 
WASP-3 was produced jTripathi et alj|20l"ol) . containing new re- 
sults including Ka = 290.5t^ | kms'^ and A = 2>.itl\ °- 

The third and most recent Ka has been used to obtain the 
physical properties of the WASP-3 system (Table A46). I find that 
WASP-3 b has a rather la rger mass and radius than most literature 
studies, except for that of iTripathi et al.l ( |2010|) . More precise mea- 
surements of Tcft and for WASP-3 would be useful to improve 
our understanding of its physical properties. More extensive spec- 
troscopy would also be useful to pin down Ka, for which a variety 
of measurements are currently available. 



4.11 WASP-10 



WASP-10 was found to be a TEP system bv lChristian et alj j2009l) . 
and is notable for containing a fairly massive (3.2 Mjup) planet 
transiting a small (0.70 Rq) and low-mass (0.75 M©) star. The 
transit events are a generous 3.0% dee p, so photometric foU ow- 
up of this system is comparatively easy, (christian et ^ ^20091) ob- 
tained data from the 0.8 m Tenagra and 1.2 m Mercator telescopes, 
but unfortunately none of the datasets cover a full transit event. 
High-precision follow-up p hotometry of o ne com plete transit of 
WASP-10 was obtained by Ijohnson et al.l ( l2009b[) , using a novel 
orthogonal frame transfer CCD (OPTIC) to shape the point spread 
function and thus obtain a scatter of only 0.5 mmag per datapoint 
at a reasonable sampling rate. I analyse the OPTIC observations to 
obtain the photometric parameters, and the Mercator data to pro- 
vide a consistency check. 

One complication for WASP-10 is its eccentric orbit. This 
is handled in the way described in Sect. |2.3l using the con- 
straints ecostj = -0 .045 ± 0.02 and esinu = 0.023 ± 0.04 
jjohnson et l2009bl) . For both light curves I find that corre- 
lated noise is unimportant and that the LD fit/fix solutions are 
best (Fig.llSt. The two sets of results agree well (Table A49) and 
I adopt the OPTIC ones as final. The values agree well with those 
of Johnson et al. (2009b) but my errorbars are rather larger, due in 
part to the inclusion of several different LD laws rather than the 
reliance on only one. The agreement with lChristian et all ( |2009|) is 
less good but still acceptable. 

The ABSDIM analysis returns res ults (Table A50) which are 
again in good agreement with those of Ijohnson et al.l j2009bh but 
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Figure 15. Phased light curves of WASP-10 compared to the best fit found 
using JKTEBOP and the quad ratic LD law. The u pper curve and residuals 
represent the Mercator data ^Christian et al.| |2009l) and the lower curve and 
residuals are the OPTIC data jjohnsoneran[2009hl) . All offsets are additive 
in flux. 



with larger errorbars. The agreement between different model sets 
is unusually poor for WASP-10, resulting in systematic errorbars 
which are a significant fractio n of the random errorbars and as large 
as the total errorbars given by Johnson etjil.; Table A50 provides the 
first measurement of O for WASP-10 and also corrects a calcula- 
tion error in the Tjq listed in one of the published studies of this 
system. An improved photometric study of WASP-10 would settle 
the existing disagreement on its k value, and the system would also 
be favoured by more precise T^a and [|^] measurements. 

After the above study of WAS P-10 was completed, new 
datasets on thi s star w ere presented bv lOittmann et alj ( 1201 Oh and 
iKrejcova et al ]( l2010h and a disagreement over the system proper- 
ties became manifest. T he two new studies both prefer system prop- 
erties similar to those of Ichristian et alj j2009l) but d iscrepant with 
the resu lts - based on a much better light curve - of Ijohnson et al.l 
( l2009bl) . This will be revisited in the future, once the newer data 
and perhaps further observations become available. 



4.12 XO-2 



The s econd planet discovered by the XO project dMcCullough et all 
l2006h is noteworthy for having a parent star which is very metal- 
rich ([^] ~ 0.4 5) and a me mber of a common proper mo- 
tion binary system teurke et alj|2007h . Good transit light curves 
from the Perkins telescop e were published in the discovery paper 
dMcCullough et al.ll20odl and from the FLW O 1.2 m by the Tran- 
sit Light Curve Project jpemandez et alj200^ . The s ystem has also 
been observed as part of the NASA EPOXI mission dBallard et al.l 
|2009). 

In this work I analyse the Perkins and FLWO observations 
(Fig.ll6t. In the former case correlated noise is important and in 
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Figure 16. Phased light curves of XO-2 compared to the best fits found 
using JKTEBOP and the qua dratic LP law. The upper light curve is i?-band 
from the Perkins telescope I Burke et al.||2007 [) and the l ower one is z-band 
from the FLWO 1.2 m telescope liF ernandez et alj|2009l) . The residuals are 
offset from zero to the base of the figure. 



the latter it is not. I adopt the LD fit/fix results for both datasets 
and combine their probabihty density functions to find the final re- 
sults. The solutions have i ^ 90° and thus asymmetric errorbars. 
The agreement between the two light curves and versus published 
values is good (Table A5 3). 

The ABSDIM analysis is complicated by the high [^] (Ta- 
ble A54). The results using the five different theoretical models are 
scattered, giving systematic errors which are larger than the statisti- 
cal ones for the two most-affected quantities, Ma and a. However, 
the agreement with other studies is good. A more precise Ka value 
would be profitable. 



4.13 XO-3 

XO-3 was discovered by Ijohns-Krull et alj bOOSi) to be a TEP 
which is so massive it is near the 13 Mjup value which represents 
the minimum mass of a brown dwarf. Ijohns-Krull et ai]( l200a pre- 
sented two alternative sets of physical properties for the system, 
the first of which put X0-3b at 13.25 ± 0.64M,Tup but yielded 
a relatively poor fit to the observed transit light curve. The sec- 
ond set ignored the spectroscopic measurement of loggA, yield- 
ing M2 = 12.03 ± 0.46 Mjup and a much better fit to the light 
curve. The latter solution is preferable because light curve shapes 
are more reliable than spectroscopically-derived surface gravities 

for late-type dwarf stars. 

A follow-up study of XO-3 was presented by Iwinn et all 

( l2008bh , based on high-quality new light curves and also ignoring 
the spect roscopic log gA. Additional observations have been pre- 
sented by Winn et al. ( 20 0j[). In th is work I analyse the 2-band pho- 
tometry obtained by Winn et al. I l2008b) using the FLWO 1.2 m, 
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Figure 17. Phased light curves of XO-3 compared to the best fits found 
using JKTEBOP and the quadr atic LD law . From top to bottom the light 
curves are KeplerC am ^-band jWin ri et al.ll2Q 08b). and KeplerCam r-band 
and FLWO /-band jwinn et al.l2009l) . The residuals are offset from ze o. 



and t he r-band FLWO 1.2 m and J-band Nickel datasets presented 
by Winn et alj ll2009h. The photometric observations presented by 
Ijohns-KruU et al.l l l2008h are not included because they comprise 
many datasets of only moderate quality. 

An important property of XO-3 is its substantial orbital ec- 
centricity (e = 0.28), which is a common feature of the more mas- 
sive planets jSouthworth et al. 2009c) and might indicate that they 
are a different population of objects to their less massive cousins. 
XO-3 is also known to exhibit a larg e spin-orbit misalignment 
l lHebrardet^ l2008l : IWinn et alj 120091) sugges tive of dynamical 
evolution through gravitational interactions (e.g. lFabrvckv & Wintil 
2009). In the following analysis I adopt the const raints e = 
0.2884 ± 0.0035 and ui = 346.3 ± 1.3, taken from IWinn eFal] 
(i2009i) . Compared to a solution assuming a circular orbit, ta and 
Tb decrease by 0.0033 and 0.00029 respectively, which is less than 
la in both cases. 

Correlated errors are unimportant for all three light curves, 
and in each case the best solutions are LD fit/fix. The z and I data 
agree well but the r results have a higher i and a 2.5a lower t-a 
(Table A58). The most extensive dataset is z, so I have combined 
the results from this and I and rejected the r results as discrepant. 
T he best fits are plott ed in Fig|17| My results agree well with those 
of lWinnetal] ( l2009h and with the second of the two alternative 
solutions given bv ljohns-Krull et al.l ( l2008l) . 

When determining the physical properties of the XO-3 sys- 
tem I adopted increased uncertainties of 75 K in Tcff and 0.05 dex 
in , to allow for the possibility of systematic errors in these 
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Figure 18. Phased ligh t curve of XO-4 from the Perkins telescope 
iMcCullough et al.ll20O8h . The blue line shows the best fit from JKTEBOP 
using the quadratic LD law. The residuals are offset from zero. 



values for low-mass stars such as XO-3 A, and to account for the 
low spectroscopic gravity valu e (3.95 ± .06 versus 4.23 ± 0.04) 
in the discovery paper ( Johns -Krull et al . 2008). Like the JKTE- 
BOP out come, the re sults of the ABSDIM analysis agree well with 
those of IWinn et alj (i2009i) but not with the preferred solution of 
Ijohns-Krull et al.l feOOsir The mass of the planet is A/b = 11.8 ± 
0.5 Mjup, which is close to but below the 13 M,jup dividing line be- 
tween planets and brown dwarfs. The VRSS model results disagree 
strongly with those of the other models, so were not included in the 
final analysis. The T^q of XO-3 b is high at 1729 ± 34 K, making it 
an interesting object for the study of planetary atmospheres. Aside 
from the VRSS models, XO-3 is one of the best-measured TEPs. An 
improved spectroscopic study, incorporating the best log qa value 
given in Table A59, would be the best way of improving this under- 
standing even further. 



4.14 XO-4 



XO-4 was discovered to be a TEP by iMcCullough et al.1 ( l2008h : 
the parent star is one of the hottest of the known planetary hosts. 
H ere I analyse the J^-band Perkins telescope light curve obtained 
bylM cCullough et alJJioOSl) . Correlated errors are unimportant and 
the best solutions are LD fit/fix. The inclination is near 90° , result- 
ing in asymmetric errorbars. My photometric sokition (Fig.llSt is 
in excellent agreement with that of iMcCullough et al.l (Table A61). 

The results of the ABSDIM analysis are given in Table A62, 
and include the first reported measurements of the Tpq, Safronov 
number, gb, Pa and pb of the XO-4 system. The other output pa- 
rameters agree well with those of McCuUough et al. (200JP; Mb 
is 0.2 Mjup smaller in my solution but this is within the errorbars. 
Improved photometric and RV observations of XO-4 would be ben- 
eficial. 



4.15 XO-5 

The discov e ry tha t XO-5 is a TEP system was presented by 
iBurke et alj i2008t) and extensive follow-up observations in the 
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Figure 19. Phased light curves of XO-5 compared to the best fits found 
using JKTEBOP and the quadratic LD law. From top to bottom the light 
curves are Perkins i?-b and iBurke et al,,2Q0 8). and KeplerCam i-band and 
2 -band jPal et alj|2009l) . The residuals are offset from zero. 



context of the HAT consort ium were presented bv lPal et aL 
I adopt the Ka value from lPal et al.l ( l2p09|). which agrees with but 
is much more precise than that given bv lBurke et al. com- 



parison of the Tcff and given by the two studies shows a con- 
ceming d isagreement, which may be due to the treatment of log qa- 
IPaletalJ ( 12009,) fix log qa at the transit-derived value during their 
spectral synthesis analysis, and also have spectra with a higher sig- 
nal to noise ratio, so I have preferred their spectroscopic results. 

Three high-quality light curves of XO-5 are availa ble: J?-band 
cover age of one full transit using the Perkins telescope jBurke et al.1 
|2008|) . and i- and z-band observations (both covering one full 
and one partial transit) obtained with the FLWO 1.2 m telescope 
jPal et alj2009h . 1 have analysed all three datasets, and the best fits 
are shown in Fig.[T9] In the case of the R and z data correlated 
noise is important, and in all cases the best solutions are LD fit/fix. 
The resulting photometric parameters are in good agreement except 
for k, for which there is a large scatter of 5a (Table A66). The three 
parameter sets have therefore been combined and the errorbar in k 
increased to account for this discrepancy. 

The physical properties of XO-5 are given in Table A67, and 
are in good agreement with established values. This system is 
a good candidate for improved photometric observations, which 
would allow to sort out the discrepancy in k and improve the preci- 
sion of the system properties. 
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4.16 TEPs studied in previous papers 

The preceding subsections liave presented full studies of fifteen 
TEPs. In this subsection I apply my improved ABSDIM analysis to 
the other fifteen systems for which photom etric results have already 
been calculated tPaperl JSouthworth et al . 2009a b e, 2010). Their 
physical properties are collected and compared to literature results 
in Tables A68 to A81, and relevant points are discussed below. 

GJ 436. Attempts to obtain solutions for Tcft = 3350 K failed. 
Inspection of Fig.|4]explains this: for a 0.5 star we expect Tcft 
values in the range 3500-4000 K (except for the Claret models 
which prefe r 3100 K). I have t herefore adopted the higher T^s of 
3500 K from lBean eFZI 12006), with an errorbar doubled to 100 K. 
The results with the Claret models are discrepant so were not incor- 
porated into the final solution. This does not mean that the Claret 
models are wrong - they are in fact closer to the measured prop- 
erties of eclipsing binaries (primarily CM Dra) in this mass regime 
- but that they depart from the consensus established by the other 
model sets. The lesson here is that we require an improved under- 
standing of low-mass stars to better measure the physical properties 
of the important GJ 436 system. The revised results for GJ 436 (Ta- 
ble A68) are slightly smaller than those found in Paper II, and are 
in good agreement with literature studies. 

H AT-P-1. I used a new Ka value from Ijohnson et al.l 
( l2008l) . My ABSDIM solution s prefer the Tea value used 
in Paper II teakos et alj 12007 al) to the higher one found by 
lAmmler-vonEiffetalJj2009l) . 

H P 189733. A new Ka value is a vailable fromlBoisse et al.l 
( l2009l) . supported by the value given bv lTriaud et alj ( |2009|) . The 
ABSDIM solutions with different models converged on either a very 
young (~1 Gyr) or old (9-13 Gyr) age, resulting in large system- 
atic errorbars for most ou tput parameters . HP 189733 is a n active 
star teouchv et al.ll2005bl) with starspots jPont et al.ll2007l), a short 
rotati on period of 12 d jWinn et al.]|2007bh . RV jitter jSoisse et all 
l2009l). and CaH&K emission modulated on its rotation period 
I Moutou et alj|2007|). These facts im plv a young age cSkumanich 
1 19721 : iMamaiek & HillenbrandfioOSl) so the search for the best so- 
lution was restricted to ages below 5 Gyr. This resulted in much 
more consistent solutions, with ages of 0. 1-2.9 Gyr, which were 
accepted as the final results (Table A70). 

HD 209458. I used the revised value of Ka = 84.67 ± 0.70 
given by TWH08. The revised results are in good agreement with 
those from Pap er II. Ra has been measured directly by interfer- 
ometric means l lvan Belle & von Braunll2009h to be 2(7 larger than 
found here and in all other studies of this object given in Table A7 1 . 

OGLE-TR-10. The light curve solution used in Paper II was 
unintentionally a preliminary rather than the final one from Paper I. 
This has been corrected here, resulting in a less dense planet. The 
properties of OGLE-TR-10 are rather uncertain and the system is 
badly in need of improved spectroscopy and photometry. 

OGLE-TR-132 . I use the more pr ecise Ka value of 167 ± 
18kms~ given bv lMoutou eT^ j2004l) which was overlooked in 
Papers I and II. OGLE-TR-132b has a high T^'q of 2017 ± 97 K 
(Table 75). 

WASP-l. Th e results for this system a re reproduced from the 
dedicated study of lSouthworth et al.[|201oh . 

WASP-4, WASPS, WASP-1 8 I use the new and improved Ka 
values for these three TEPs from lTriaud et al.l j2010l) . 



5 TRACKING THE SYSTEMATIC ERRORS IN THE 
PROPERTIES OF TRANSITING SYSTEMS 

A major result of the current work is a detailed understanding of 
where model-dependent systematic errors surface in the analysis 
of TEPs, and the importance of these systematics for the various 
physical properties which can be calculated. The approach used in 
this work means that all of the model dependence is combined into 
i^b (the velocity amplitude of the planet), making it an excellent 
tracer of systematic errors. For each TEP a separate value of K\y 
is found using each of the five sets of stellar models, as well as 
for the empirical mass-radius relation (Paper II). I have converted 
these into 'consensus values', (Kh), using the same algorithm as 
for the other measured physical properties: the unweighted mean 
of the values from the five different stellar model sets. Remember 
that (Kh) should not be taken as an indicator of correctness, only 
of concordance. 

From Eqs. 4 to 1 2 in Paper II it can be seen that the component 
masses are most sensitive to systematic errors (Ma oc K^ and 
A/b oc Ky^) and that their radii, gravities and densities are less so 
(all directly proportional to Kh). The semimajor axis is also rather 
model-dependent as the other quantities required to calculate it all 
have much smaller uncertainties than A'b does. 

A detailed exploration of model-dependence is presented in 
Fig.|20]as a function of Tcs, and Ma. The top panels in this 
Figure highlight the generally good agreement between different 
model sets. The mass-radius relation is in much poorer agreement 
and is biased to high values as it does not account for the effects of 
evolution through the main-sequence phase (particularly apparent 
for HAT-P-2, TrES-4 and XO-3). This bias to large A'b comes from 
trying to reproduce the low density of an evolved star, and pushes 
the mass and radius of the planet and star to high and incorrect 
values. 

Turning to the stellar model sets in Fig.|20l clear trends with 
respect to TcB, [^] and Ma can be seen in many cases. The Claret 
models yielded larger Kh values on average, particularly for low 
values of TcS and [|^] . This predilection for high values is bal- 
anced by other model sets: the Teramo and VRSS models both tend 
to produce lower Kh values and the DSEP models trend to low A'b 
for cooler and more metal-poor stars. Within this melange, the 
models are the closest to the consensus value and do not exhibit sig- 
nificant trends with the stellar properties. The models are thus 
the best choice to obtain quick results, although of course other 
model sets are required for the assessment of systematic errors. 

Figs.l21landl22lcompare the sizes of the random and system- 
atic errors in stellar mass, as a function of Ma and of [^] , re- 
spectively. A/a was chosen for this comparison because it is one 
of the parameters more affected by systematics; it was also ex- 
pected that these systematic effects would be minimised in the re- 
gion of 1 Mq and [^] = as all of the stellar model sets are 
calibrated on the Sun. Surprisingly, this does not turn out to be the 
case. There is a hard lower limit of 1% on the errors in Ma arising 
from model-dependent systematics, but this lower limit is reached 
at a wide range of masses and [^] values. The systematic errors 
are clearly larger for lower-mass stars (<O.9M0) and for a high 
metal abundance ([§] > 0.4). The hard lower limit is unavoid- 
able until stellar model sets are in much better agreement, and the 
mass-radius relation results hint that the real systematic errors are 
probably somewhat higher than this. The random errors in Ma tend 
to decrease towards higher masses (temporarily ignoring the much 
fainter OGLE systems), as expected because more massive stars are 
intrinsically brighter and therefore easier to obtain good data for. 
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Figure 20. Comparisons between the K]^ values obtained using specific sets of stellar evolutionary models and the unweighted mean value, (K^,), for each 
TEP. From left to right the panels show results for the mass-radius relation and then the five stellar model sets. The top panels compare K^, to (K^,) for each 
model set, with parity indicated by a dotted line. Lower panels show the difference (K^, — (i^b)) as functions of effective temperature, metal abundance and 
stellar mass. Errorbars have been ignored for clarity; their median values are ±3kms^^ (statistical) and ±1.4kms^^ (systematic). 



5.1 An external test: WASP-11 versus HAT-P-10 

One of the best ways to investigate the presence of systematic errors 
is to compare two independent studies of the same object. Whilst 
multiple discovery and characterisation publications have been pre- 
sented for several TEPs (e.g. XO-1 and XO-5), successive papers 
on these objects have in each case been informed by the initial dis- 
covery papers. There are two exceptions: HD 80606 and WASP- 
11 /HAT-P-10. In the case of HD 80606 three groups discovered its 
transiting nature essentially simultaneously, but their analyses were 
heavily dependent on the same published spectroscopic observa- 
tions. WASP-11 /HAT-P-10 is therefore the only TEP which was 
discovered and fully characterised by two groups working with- 
out knowledge of each others' analyses. The discovery papers were 
submitted within a week of each other, but the agreed name of the 
system is WASP-1 1 /HAT-P-10 because the WASP group submit- 
ted their paper first. 

In Table|4]l collect t he properties of WASP-11 determin ed in 
the two discovery papers (West et al. 2009:'Ba kos et alj20 09V The 
values in general show a gratifying agreement, but two quantities 



stand out as being discrepant to some extent. The Tcff values dis- 
agree by 1.6(T and the k values by 2.7a. The divergent k values 
have a large knock-on effect on the planetary properties (Rb, gh 
and /9b). 

The possibility of systematic errors in Tcff is a well-known 
phenomenon (see Sect.[3]l, so the 1.6(t disagreement is unsurpris- 
ing. Similarly, k measurements primarily depend on the observed 
transit depth (see Sect.|2ll and can be affected by imperfect normali- 
sation as well as by astrophysical effects such as starspots. Discrep- 
ancies in k have previously been found in the we ll-studied systems 
HD 189733 and HD 209458 (Paper I), WASP-4 ( Sou thworth etahl 
l2009bi) . HAT-P-1 and WASP-10 (SectionsglJandgTI). 

This implies that a high-quality study of a TEP carmot rely on 
photometric coverage of only one transit, irrespective of its quality, 
but must include observations of two or preferably more in order 
to pick up on well-camouflaged systematic errors affecting light 
curves. It is important to realise that a systematic over- or under- 
estimation of a transit depth has a big effect on measurements of the 
planet properties, but cannot be identified in any dataset covering 
only one transit. 
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Table 5. Physical properties of the stellai' components of the TEPs studied in this work. For each quantity the first uncertainty is derived from a propagation of 
all observational errors and the second uncertainty is an estimate of the systematic eiTors arising from the dependence on stellar theory. 



System 


Semimajor axis (AU) 


Mass (Mq) 


Radius (Rq) 


log QA [cm/s] 


Density (p0) 


Age (Gyr) 



GJ436 

HAT-P-1 

HAT-P-2 

HD 149026 

HD 189733 

HD 209458 

OGLE-TR-10 

OGLE-TR-56 

OGLE-TR-111 

OGLE-TR-113 

OGLE-TR-132 

OGLE-TR-I82 

OGLE-TR-211 

0GLE-TR-L9 

TrES-1 

TrES-2 

TrES-3 

TrES-4 

WASP-1 

WASP-2 

WASP-3 

WASP-4 

WASP-5 

WASP- 10 

WASP- 18 

XO-1 

XO-2 

XO-3 

XO-4 

XO-5 



0.02887 ± 0.00089 
0.05535 ± 0.00057 
0.06740 ± 0.00074 

n nA9SS +0-00048 
0.042SS _o.oo027 

0.03142 ±0.00038 
0.04747 ± 0.00046 
0.04516 ±0.00099 
0.02386 ± 0.00028 
0.04651 ±0.00099 
0.02278 ± 0.00047 
0.03029 ±0.00062 
0.05205 ±0.00057 



0.05105 



-1-0.00076 



■0.00073 

0.0404 ±0.0011 

0.03946 ± 0.00039 

0.03635 ± 0.00063 
n n99S'? +0.00012 

U.U^^BJ -0.00025 

0.04965 ± 0.00087 



0.03898 



-1-0.00036 



0.00039 

0.03033 ±0.00060 
0.03187 ±0.00086 

+0.00053 
0.00055 

0.02714 ±0.00049 

0.0378 ±0.0014 

0.02034 ±0.00032 

0.04944 ±0.00062 

03647 +0 00059 
U.UJO^/ -0.00058 

0.04529 ±0.00057 

0^5475 +0 00094 
U.U04I-0 _o. 00051 

0.0494 ±0.0011 



0.02307 ' 



± 0.00035 
± 0.00042 
±0.00034 

+0.00013 
-0.00019 

±0.00036 
±0.00031 
±0.00015 
± 0.00009 
±0.00051 
± 0.00022 
±0.00018 
±0.00031 

+0.00020 
-0.00021 
± 0.0002 
± 0.00060 

± 0.00035 

+0.00012 
-0.00017 

± 0.00028 

+0.00013 
-0.00014 

±0.00043 
± 0.00020 

+0.00026 
-0.00022 

± 0.00022 
±0.0008 
±0.00016 
± 0.00050 

+0.00069 
-0.00081 

±0.00045 

+0.00022 
-0.00035 

±0.0002 



0.459 ± 0.043 
1.134 ±0.035 
1.279 ±0.042 

, 97-1 +0.043 
^•^'^ -0.024 

0.840 ±0.030 
1.148 ±0.033 
1.277 ±0.082 
1.233 ±0.043 
0.833 ±0.054 
0.768 ± 0.048 
1.297 ±0.078 
1.187 ±0.039 

-. q-. r, +0.059 
i.oiz _o.o56 

1.42 ±0.11 
0.892 ±0.026 
1.049 ±0.054 



±0.017 
± 0.026 
±0.019 

+0.012 
-0.017 

± 0.029 
± 0.022 
±0.013 
±0.014 
± 0.027 
± 0.022 
±0.023 
±0.021 

+0.015 
-0.016 
±0.02 



0.454 ± 0.029 ± 0.005 
1.1 12 ±0.031 ±0.008 
1.68 ±0.15 ±0.01 

1 9Qn +0-120 +0.004 
l.ZMU _o.o58 -0.006 

0.752 ± 0.023 ± 0.009 
1.162 ±0.012 ±0.008 

1.52 ±0.10 ±0.00 
1.26 ±0.14 ±0.00 
0.842 ± 0.042 ± 0.009 
0.780 ± 0.029 ± 0.008 
1.37 ±0.14 ±0.01 

1.53 ±0.17 ±0.01 



1.56 



+0.18 +0.01 



090+0-014 
u.aza _o.030 



± 0.030 

+0.014 
0.021 

1.292 ±0.067 ±0.022 

+0.013 
0.014 

± 0.034 



0.10 -0.01 

1.503 ±0.083 ±0.008 
±0.041 0.818 ±0.017 ±0.013 
1.002 ±0.029 ±0.010 



n SIS +0-011 +0.004 

u.oio _o.oi3 -0.006 

1.92 ±0.11 ±0.01 



1 94-? +0-034 - 
l.z^o _o.o37 

0.803 ± 0.049 
1.26 ±0.10 
0.914 



1.455 



+0.052 +0.005 



-0.064 
-0.066 

1.004 ±0.055 
0.752 ±0.081 
1.256 ±0.059 
1.037 ±0.039 
946 

-0.046 

1.206 ±0.046 



±0.02 

+0.031 
-0.026 
± 0.025 
± 0.048 
± 0.029 



±0.031 

+0.054 
-0.062 

± 0.036 

+0.016 
0.024 

0.914 ±0.064 ±0.010 



-0.079 -0.005 

0.807 ±0.019 ±0.011 

1.377 ±0.085 ±0.009 
n ont^ +0-021 +0.010 

u.auo _o.022 -0.009 
1.077 ±0.042 ±0.009 
0.703 ±0.036 ±0.015 
1.222 ±0.042 ±0.010 
0.942 ± 0.022 ± 0.009 



4.787 ± 0.030 
4.400 ± 0.024 

4.092 ± 0.074 

4 321+0-042 

4.610 ±0.026 
4.368 ±0.005 
4. 178 ±0.053 
4.331 ±0.094 
4.508 ± 0.044 
4.539 ±0.028 
4.275 ± 0.083 
4. 142 ±0.089 



4.170 



+0.052 



1 9S"^ +0-068 
i.ZSO _o.036 



n Q7n +0-028 +0.018 
U.y /U _o 035 -0.022 

1.409 ±0.054 ±0.014 

1 ^-^n +0-362 +0.006 
l.oou _o.o69 -0.010 



-0.085 

4.236 ±0.043 
4.563 ±0.019 
4.457 ± 0.027 

A t^SI +0-007 
4.581 „o.010 

3.981 ±0.047 

4 207+0-045 
4.ZUI _o.o28 

4.529 ±0.017 
4.262 ± 0.044 

A /tat; +0.011 

4.485 _o 012 
4.375 ±0.030 
4.620 ± 0.049 
4.363 ±0.027 

4.506 ±0.021 

4 440+0-037 

„o.021 

4.222 ± 0.027 



±0.005 
± 0.003 

± 0.002 

+0.001 
-0.002 
±0.005 

± 0.003 
±0.001 
± 0.002 
± 0.005 
±0.004 
±0.003 
±0.003 

+0.002 
-0.002 

±0.002 
±0.007 
±0.004 

+0.002 
-0.003 
± 0.002 
+0.002 
-0.002 

± 0.006 
±0.003 

+0.005 
-0.004 

± 0.004 
± 0.009 
±0.003 



± 0.004 

+0.008 
-0.010 

±0.004 

+0.002 
0.002 

1.065 ± 0.064 ± 0.004 4.344 ± 0.043 ± 0.002 



4 170+0.034 - 



4.92 ± 0.55 

0.824 ± 0.066 

0.268 ± 0.070 
592+0-083 

1.98 ± 0.17 
0.733 ± 0.008 
0.361 ± 0.063 
0.62 ± 0.21 
1.40 ±0.19 
1.62 ± 0.13 
0.50 ±0.15 

0. 33 ± 0.10 

345+0-068 
u.oij_o 090 

0.418 ±0.061 
1.632 ± 0.092 
1.043 ±0.088 

1 700+°-''*^ 

1. /UU_o.o51 

0.182 ± 0.030 

403+0-069 
u.'iuo_o 037 

1.527 ±0.067 
0.484 ± 0.073 

, r,oq + 0.020 

i.zJO_Q 022 

0.803 ± 0.080 

2.16 ±0.31 
0.689 ± 0.062 
1.242 ± 0.080 

1 n'^7+0-1'^^ 

l.uoi _o.058 

0.431 ± 0.041 
35Q+0-046 

0.76 ±0.11 
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Figure 21. Plot of the sizes of the statistical (red filled circles) and system- 
atic (blue crosses) errorbars for the TEPs studied in this work, versus the 
stellar masses. The errors are plotted as fractions of and the statistical 
and systematic errors for each TEP are connected by grey lines. 



6 PHYSICAL PROPERTIES OF THE TRANSITING 
EXTRASOLAR PLANETARY SYSTEMS 

The main result of this work is the determination of the basic physi- 
cal properties of thirty TEPs using homogeneous methods. Detailed 
results for each TEP are given in the many Tables in the Appendix, 
and the final results for all systems have been summarised in Tables 
[5]and[6l The homogeneous nature of these numbers means they are 
well suited for comparing different TEPs, for planning follow-up 
observations, and for performing detailed statistical studies. 

The masses and radii of the stars considered in this work are 
plotted in Fig.[23]with both their random and systematic errorbars, 
and with the empirical mass-radius relation (Paper II) indicated. A 
striking aspect of Fig.[23]is that for Ma < 1.2 M© the stars gen- 
erally stick close to the main sequence, but beyond this point their 
distribution turns almost vertically upwards. The two exceptions 
are XO-5, which has a more evolved lower-mass star, and OGLE- 
TR-L9, which contains a more massive star (1.4 M©). 

That low-mass stars stick closely to the main sequence is to 
be expected due to their long evolutionary timescales. The shift to 
larger Ra at 1.2 M0 is reasonable because a larger stellar radius 
makes it more likely that a given planet is transiting (see Paper II). 
The avoidance of more massive stars is likely due to difficulty of 
measuring RVs for these objects, making them lower-priority and 
more observationally expensive targets for TEP search consortia at 
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Table 6. Physical properties of the planetaiy components of the TEPs studied in this work. For each quantity the first uncertainty is derived from a propagation 
of all observational eiTors and the second uncertainty is an estimate of the systematic errors arising from the dependence on stellar theory. 



System 



Mass(Mjup) 



Radius (Rjup) 



9b (ms ^) 



Density ( pjup) 



To'q (K) 







GJ436 

HAT-P-1 

HAT-P-2 

HD 149026 

HD 189733 

HD 20945 8 

OGLE-TR-10 

OGLE-TR-56 

OGLE-TR-Ul 

OGLE-TR-113 

OGLE-TR-132 

OGLE-TR-182 

OGLE-TR-211 

0GLE-TR-L9 

TrES-1 

TrES-2 

TrES-3 

TrES-4 

WASP-1 

WASP-2 

WASP-3 

WASP-4 

WASP-5 

WASP- 10 

WASP- 18 

XO-1 

XO-2 

XO-3 

XO-4 

XO-5 



0.0737 ±0.0051 ±0.0018 
0.524 ±0.016 ±0.008 
8.74 ±0.25 ±0.09 

n +0.013 +0.002 
u.oao _o.oil -0.003 



1.150 

0.714 

0.68 

1.300 

0.54 

1.24 

1.17 

1.06 



± 0.028 

±0.014 

±0.15 

± 0.080 

±0.10 

±0.17 

±0.14 

±0.15 



± 0.027 

± 0.009 

±0.00 

±0.010 

±0.01 

±0.02 

±0.01 

±0.01 



p. rjr- +0.15 +0.01 
''^ -0.15 -0.01 

4.34 ±1.47 ±0.05 
0.761 ±0.045 ±0.023 
1.253 ± 0.047 ±0.024 

1 Q1 n +0.060 +0.020 
I.MIU _Q 070 -0.029 

0.877 ±0.072 ±0.010 
n sen +0-072 +0.006 

u.ouu _Q 072 -0.006 

0.847 ±0.038 ±0.024 
2.06 ±0.13 ±0.03 

1 907 +0.059 +0.028 
-0.062 -0.023 

I. 565 ±0.058 ±0.026 
3.16 ±0.23 ±0.13 
10.29 ±0.33 ±0.16 
0.924 ±0.075 ±0.019 

n t;>^^ +0.062 +0.021 

u.ooo _o,o57 -0.025 

II. 83 ±0.31 ±0.23 

1 c;91 +0.160 +0.012 
l.OZi _Q _o.019 

1.084 ±0.055 ±0.008 



0.365 ±0.018 
1.217 ±0.038 

1.19 ±0.12 

610+0 099 
1.151 ±0.036 
1.380 ±0.015 
1.72 ±0.11 

1.20 ±0.17 
1.077 ±0.072 
1.111 ±0.049 
1.25 ±0.16 
1.47 ±0.14 

1 9fi9 +0' 158 
l.zuz _o.o9l 



± 0.004 
± 0.009 

±0.01 

+0.002 

-0.003 

±0.013 

± 0.009 

±0.01 

±0.00 

±0.012 

±0.011 

±0.01 

±0.01 

+0.005 
-0.005 

± 0.009 



1.614 ±0.083 
1.099 ±0.031 ±0.017 
1.261 ±0.039 



l.OUO _o.o25 

1.81 ±0.15 

1 484+0-059 
l.^O^i _o.o91 

1.043 ±0.029 
1.454 ±0.083 

, oc;7 +0.033 
1.001 _o.033 



±0.012 

+0.007 
-0.010 

±0.01 

+0.005 
-0.006 

±0.015 
± 0.009 

+0.015 
-0.013 

±0.009 



1.164 ±0.056 
1.067 ±0.064 ±0.022 
1.158 ±0.054 
1.206 ±0.039 



002 +0-034 

1.248 ±0.047 

1 00+0-38 
^■^^ -0.06 

1.089 ±0.082 



± 0.009 
±0.012 

+0.019 
-0.022 

±0.012 

+0.01 
-0.01 

±0.004 



13.7 ± 1.1 
8.77 ±0.56 
152 ± 30 

9Q 7+6.8 
^•J- ' -6.2 

21.5 ± 1.2 

9.30 ±0.08 

5.7± 1.4 

22.3 ± 6.7 

11.5 ± 2.5 

25.0 ±3.7 

18.5 ± 5.0 

12.1 ± 2.9 

ii-6l^;^3 

41 ± 14 

15.6 ± 1.2 
19.5 ± 1.1 

6.7± 1.2 

19.32 ± 0.80 

24.2 ± 2.6 
16.65l»;f3 

28.7 ±2.6 
68.9 ± 6.7 

190 ± 16 

15.8 ± 1.5 
14.01?- 
188 ± 13 
22.813-2 

22.7 ±3.2 



1.51 ±0.19 
0.290 ± 0.027 
5.1 ±1.5 



1.57 



fO.72 
-0.57 

0.755 ±0.066 
0.272 ± 0.004 
0.134 ±0.038 
0.75 ±0.34 
0.43 ±0.11 
0.91 ±0.16 
0.59 ±0.24 
0.332 ±0.1 11 

U.OIZ _o,i32 
1.03 ±0.37 
0.573 ±0.056 
0.625 ±0.051 

0.8601^^7 
0.148 ±0.039 

263+°-''5** 
u.zuo _o.036 

0.747 ± 0.047 

0.67 ±0.11 



0.495 



+0.015 



-0.018 

0.99 ±0.14 

2.60 ±0.39 

6.64 ±0.90 

0.526 ±0.063 

569 +"-"9 
u.aua _o.o72 

6.08 ±0.67 

71+0-13 
'J-'l -0.40 

0.84 ±0.18 



±0.02 
±0.002 

±0.0 
+0.01 
-0.01 
±0.009 
± 0.002 
±0.000 
±0.34 
±0.00 
±0.01 
±0.00 

±0.002 
+0.002 
-0.001 
±0.01 

±0.009 

±0.006 

+0.007 
-0.004 

±0.001 

0.001 
0.001 

±0.010 
±0.00 

+0.005 
-0.006 
±0.01 
±0.05 
±0.05 
±0.005 

+0.013 
-0.011 

±0.06 

+0.01 
-0.00 
±0.00 



669 ± 22 
1291 ±20 
1516 ±66 

1626l|9 

1191 ±20 
1459 ± 12 
1702 ± 54 
2140 ± 120 
1034 ± 28 
1355 ± 35 
2017 ±97 
1550 ± 81 
16861^° 
2039 ± 51 
1147 ± 15 
1467 ± 27 

1630122 
1861 ± 54 
18001^2 

1281 ± 21 
2028 ± 59 

I661I30 
1732 ±41 

972 ± 31 
2392 ± 51 
1210 ± 16 

13281^^ 
1729 ± 34 
1630lif 
1203 ± 33 



0.0253 ±0.0015 
0.0419 ±0.0016 
0.771 ±0.077 
0393+0-0054 

U.U030 -0,0056 

0.0747 ± 0.0024 
0.0427 ± 0.0005 
0.0279 ± 0.0062 
0.0418 ±0.0065 
0.056 ±0.011 
0.0664 ± 0.0090 
0.0436 ± 0.0072 
0.0628 ±0.0109 
0460+°-''"97 

U.U-iDU _o,0104 



± 0.0003 
± 0.0003 

± 0.004 

+0.0002 
-0.0001 

± 0.0009 
± 0.0003 
±0.0001 
± 0.0002 
±0.001 
± 0.0007 
± 0.0003 

± 0.0004 
+0.0002 
-0.0002 
±0.001 



0.153 ±0.052 
0.0612 ±0.0037 ±0.0009 
0.0688 ± 0.0024 



0.0719 _o.oo26 
0.0373 ± 0.0042 

Q 03(50 +0.0038 
u.uouo _o.0033 

0.0613 ± 0.0021 

0.0715 ±0.0048 

0460+°-°"l2 

U.UtDU _o,0013 

0.0726 ± 0.0035 
0.298 ±0.019 
0.288 ±0.014 
0.0730 ± 0.0062 

0432 +0-0049 
U.U+OZ _o,o045 

0.711 ±0.027 

0.1006 IH^^^ 
0.1075 ±0.0082 



± 0.0007 

+0.0026 +0.0006 
-0.0004 

±0.0002 

+0.0001 
-0.0001 

± 0.0009 



± 0.0004 

+0.0004 
-0.0005 

±0.0006 
± 0.006 
± 0.002 

± 0.0007 

+0.0010 

-0.0008 

±0.007 

+0.0006 
-0.0004 

± 0.0004 



Table 4. Physical properties of the WASP-1 1 /HAT-P-1 4 system deter- 
mined by the two discovery papers. Aside from the orbital period, quantities 
without uncertainties were calculated from the results given in the respec- 
tive papers. 



Parameter 



Westetah (2009) Bakos et al. ^2009) 



PoA (d) 
Teff (K) 
Kiim (kms^i) 

mm (dex) 

i^A(ms-l) 



a(AU) 
rA 

k 

in 

Ma(M0) 

log 3A [cgs] 
Mb(Mjup) 

(Rjup) 

Sb (ms-2) 

Ph (Pjup) 
Te'q (K) 



3.722465 
4800 ±100 
<6 
0.0 ±0.2 
82.1 ±7.4 
fixed 
0.043 ±0.002 
0.0801 
0.01011 
0.1273 
89.8 



+0.0011 

-0.0008 
+0.2 
0.8 



0.77 



+0.10 



0.74'' 



-0.08 
1+0.04 
-0.03 

4.45 ±0.2 

0.53 ±0.07 
91+0-06 

'-'-^^-0.03 
14 5+1-'' 
3+0.07 
-0.11 

960 ± 70 



3.7224747 
4980 ±60 
0.5 ±0.2 
0.13 ±0.08 
74.5 ±1.8 
fixed 
0.0435 ± 0.0006 
0.0842 ±0.0019 

0.01104 
0.1313 ±0.0010 
88.6j 
0.83 ±0.03 

0. 79 ±0.02 
4.56 ±0.02 

0.460 ±0.028 

1. UUO_o 027 

12.0 ±0.8 
0.479 ± 0.042 
1020 ± 17 
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-0.4 
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Figure 22. Same as Fig. l21l except that the error sizes are plotted versus 
stellar rather than mass. 



the follow-up stage. There may be a large and presently neglected 
population of TEPs around more massive stara^H 



11 Whilst RV surveys mostly concentrate on F, G and K stars, they have 
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Figure 23. Plot of the masses versus the radii of the stars in the thirty TEPs 
studied in this work. The statistical uncertainties are shown by black open 
diamonds and the systematic uncertainties by red filled diamonds. The em- 
pirical mass-radius relation from Paper I is shown with a blue line. 
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Figure 25. Mass-radius plot for the known transiting extrasolar planets. 
Those objects studied in this work are shown with (blue) filled circles and 
numbers taken from the literature with (red) crosses. The four gas giant 
planets in our Solar System are denoted with (black) filled circles. Dotted 
lines show loci where density equals 1.0, 0.5 and 0.25 pjup. 




Figure 24. Plot of the masses versus the radii of the planets in the thirty 
TEPs studied in this work. The statistical uncertainties are shown by black 
open diamonds and the systematic uncertainties by red filled diamonds. 
Blue dotted lines show where density is 1.0,0.5 and 0.25 p jup . 



Another notable feature of Fig.[23]is that the systematic errors 
(red filled diamonds in the Figure) are not negligible compared to 
the random errors, and in some cases (HD 189733, TrES-1, TrES-3, 
XO-2) are of a similar size or even larger. 

A similar diagram for the planets (Fig.l24t shows a much 
larger scatter in the planet properties, but also that systematic er- 



rors are much less important for these objects. The two outliers 
with very low densities are OGLE-TR-10 and TrES-4. 



6.1 Physical properties of all Isnown TEPs 

The physical properties of the thirty TEPs studied in this work have 
been augmented with the literature values for the other 41 known 
systems (as of 31/05/2010) to yield a larger but inhomogeneous 
sample. Fig.|25l shows their masses and radii on logarithmic axes, 
alongside those of the four Solar system gas giants. The dominant 
population of known TEPs continues to reside in a cloud of points 
centred roughly on l.OMjup and 1.2R,jup. There is a clear tail of 
high-mass planets, culminating in the brown dwarf CoRoT-3, as 
well as an increasing number of systems of much smaller mass and 
radius, typified by GJ 1214 and the first possible rocky exoplanet, 
CoRoT-7. 

Correlations have prev iously been notice d between Pprb an d 
gib dSouthworth et al.l2007bl) and Porb and Afb ('Maze h et al.l2005[) . 
The relevant plots are shown in Figs. |26] and [27] In both cases 
there are ten planets whose properties put them outside the range 
of Porb shown in these Figures. Neglecting this population of mas- 
sive planets (which ma y be a fundamentally different population 
to the lower mass ones: lFabrvck v & Winn' 2009: Southwor th et al] 
l2009cl) . the rank correlation test of Spearman ( 1904) returns a prob- 
ability of 99.80% (S.lcr) that the Porb-gb correlation is real and 
97.6% (2.3(t) that the Porb-^^^b correlation is real. The correspond- 
ing figures from Paper II, based on 44 TEPs, are 98.3% and 97.5%. 
Two concerns with these correlations are apparent. Firstly, they are 
insignificant if the longer-period planets are not rejecte43 Sec- 
ondly, the period-mass correlation has not increased in significance 
with the a ddition of 21 new TE Ps since Paper II. 

.Hansen & Barmai] ( l2007b divided up eighteen of the twenty 



not b een neglecti ng ones more massive than thi s; see 'Galland et all j2005l) . 
Ijohnson et al. l2007l)lLagrange et alj j2009l) and lBowler et aU 1 ,2010) . 



^■^ Merrill's theorem should be borne in mind: "When one thi'ows away 
discrepant results, the remainder are found to agree well." 
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Figure 26. Plot of the orbital periods versus the surface gravities of known 
TEPs. The objects studied in this work are shown with (blue) filled circles 
and other objects with (red) open circles. 
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Figure 28. Plot of equilibrium temperature versus Safronov number for the 
full sample of planets. Objects shown with (blue) circles were studied in this 
work; those which are just (red) errorbars were not. Th e dotted line shows 
the separation between Class A and Class B proposed bv lHansen & BarmanI 
420070. 
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Figure 27. Same as Fig. l26l but for planetary mass instead of g^,. 



TEPs then known into two classes based on their position in a di- 
agram of O versus Tcq. An updated version of the diagram can 
be seen in Fig.[28] and agrees with previous conclusions (Paper II) 
that the division between the classes was blurred into insignifi- 
cance. A dotted line at G = 0.055 has been drawn to show the ex- 
pected boundaries between Class I (0 ~ 0.07 ±0.01) and Class II 
(Q ~ 0.04 ± 0.01). There are several other equally possible lines 
which could be drawn through this diagram, which does not in- 
spire confidence in their reality. It is probable that the addition of 
new TEPs in the future will fill out this diagram with objects drawn 
from a single distribution. 



7 SUMMARY 

Measurements of the physical properties of many transiting extra- 
solar planets are available in the literature, but are not directly com- 
parable as they have been derived by different researchers using 
a variety of methods. In this series of papers I aim to derive the 
physical properties of the known TEPs using a rigorous and ho- 
mogeneous approach, resulting in quantities which are suitable for 
statistical analysis. Each TEP is tackled in two steps, the first being 
analysis of its transit light curves and the second being the determi- 
nation of its physical properties. 

The transit light curves are studied using the JKTEBOP code, 
which models binary systems using biaxial spheroids. I pay spe- 
cial attention to the limb darkening of the parent star; solutions are 
obtained for every light curve using each of five different limb dark- 
ening laws. Orbital eccentricity is accounted for by including e and 
u) as fitted parameters constrained by values and uncertainties mea- 
sured from radial velocity observations of the parent star. Parame- 
ter errors are assessed using Monte Carlo simulations, a residual- 
permutation algorithm, and the inter-agreement between solutions 
with different limb darkening laws. 

In this paper I present analyses of published transit light 
curves of fifteen TEPs (Table[2]l. Three of these (TrES-2, TrES-4 
and WASP-2) have fa int companion stars within 1.6" on the sky 
l lDaemgenet"ai]|2009l) . If neglected, this 'third light' dilutes the 
transit depth and causes systematic errors in parameters measured 
from the light curves. As a guideline, Lz ~ 5% causes the radius of 
a representative planet to be underestimated by 2%. I show that it is 
not possible to detect third light in synthetic light curves represen- 
tative of good ground-based observations. A third light value must 
instead be measured and used to constrain the light curve solutions. 

In the second stage of the analysis of each TEP, 1 augment the 
light curve results with measurements of the velocity amplitude, 
Tcft and for its parent star. I then interpolate within tabu- 
lated predictions from theoretical stellar models to find the over- 
all best-fitting physical properties of the star and the planet. This 
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Table 7. Summary of which types of additional observations would be use- 
ful for the thirty TEPs studied in this series of papers. * denotes where 
additional data would be useful, and indicates where it would be useful 
but difficult to either obtain or interpret. 



System 


Photometric 


Kauial 


Spectral 




observations 


velocities 


synthesis 


GJ436 








HAT-P-1 


★ 






HAT-P-2 








HD 149026 








HD 189733 








HD 209458 








OGLE-TR-10 








OGLE-TR-56 








OGLE-TR-111 




** 




OGLE-TR-113 


** 


** 


* 


OGLE-TR-132 








OGLE-TR-182 








OGLE-TR-21 1 




** 




OGLE-TR-L9 


* 


** 




TrES-1 


* 


★ 




TrES-2 








TrES-3 








TrES-4 


* 






WASP-1 


* 


★ 












WASP-3 








WASP-4 








WASP-5 


* 






WASP-10 


* 




* 


WASP-18 


* 






XO-1 


* 






XO-2 




★ 




XO-3 






* 


XO-4 


★ 


★ 




XO-5 






* 



model dependence causes systematic errors in the resulting quanti- 
ties, which I assess by comparing the individual results found using 
five independent sets of theoretical model predictions. This leads 
also to 'consensus values' for the physical properties; the Yonsei- 
Yale models are overall the closest to this consensus. 

Measurements of three physical properties are exempt from 
this model dependence: the planetary surface gravity can be deter- 
mined from only observed quantities; the density of the star has 
only a negligible dependence on stellar models; and the planet's 
equilibrium temperature does not depend on stellar models but does 
rest on the effective temperature scale of low-mass (F, G and K) 
stars. The model dependence specifies a minimum level of uncer- 
tainty for other physical properties, and this ranges from 1% for 
Ma to 0.6% for Mb and 0.3% for Ra, Rh and a. Reducing these 
minimum levels will require improvements in our understanding 
of low-mass stars. An external test of systematic errors was ob- 
tained by comparing the discovery papers of the WASP-1 1 / HAT- 
P-1 system. The agreement between the two studies is good for 
most parameters, but less so for the transit depth and the stellar 
Tcff • 

Tables [5] and |6] summarise the physical properties of thirty 
TEPs: the fifteen with light curve analyses in the current work 
plus another fifteen with photometric analyses in Paper I and later 
works. The statistical errors in these quantities are calculated using 
a perturbation analysis which returns complete error budgets for 



each output parameter. These error budgets indicate which type of 
follow-up observations are the best way of improving the parameter 
measurements for each TEP. In most the quality of the light curve 
is the primary limitation on measurements of the properties of the 
planet; many systems would also benefit from additional spectro- 
scopic observations. This is summarised in TableQfor the benefit 
of follow-up programmes. 

The measured properties of the 30 TEPs were then reinforced 
with those of the other 35 known objects and plotted in several pa- 
rameter diagrams. I find that the correlation between orbital period 
and planetary surface gravity is significant at the 3.1(7 level once 
clear outliers are thrown out - this significance has increased with 
the addition of new systems since Paper II. The period vs. plane- 
tary mass correlation has a significance of only 2.3(t, and this not 
increased since Paper II. Similarly, the distinction between Class I 
and Class II planets in the To'q vs. O diagram is blurring into irrel- 
evance. 

In this work I have now treated 30 TEPs, including all cur- 
rently known from the OGLE, TrES and XO surveys. The discov- 
ery rate of transiting planets is expected to continue rising, as the 
ground-based surveys HAT and WASP mature and new areas of pa- 
rameter space are probed by the space missions CoRoT and Kepler. 
A homogeneous study of the atmospheric parameters of all TEP 
host stars would be very complementary to this work. Whilst the 
light curves obtained by CoRoT and Kepler are of extremely high 
quality, it should not be forgotten that extensive velocity observa- 
tions are also needed to study transiting systems in detail. Precise 
physical properties of transiting planets are a fundamental require- 
ment of constraining the atmospheric characteristics, formation and 
evolution of planetary systems. 
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